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The  overall  objective  of  the  research  is  to  develop  and  verify 
mathematical  models  of  delamination  and  transverse  fracture  which  account 
for  local  (crack  tip)  and  global  damage  distributions.  One  specific 
objective  is  to  demonstrate  theoretically  and  experimentally  that  "work 
potentials"  (which  are  analogous  to  strain  energy)  exist  for  composites 
with  constant  and  changing  damage  and  with  viscoelastic  behavior.  The 
second  objective  is  to  develop  and  verify  methods  of  analysis  for 
predicting  crack  growth  in  elastic  and  viscoelastic  composites  with 
distributed  damage;  whenever  they  are  justified,  work  potentials  will  be 
used  to  characterize  material  behavior  in  order  to  simplify  fracture 
analysis. 

2.  STATUS  OF  THE  RESEARCH 

2.1  Introduction 

Methods  of  deformation  and  fracture  characterization  and  prediction 
are  simplified  when  strain  energy-like  potentials  based  on  mechanical  work 
can  be  t,sed,  as  described  in  the  publications  in  the  Appendix.  With  these 
so-called  work  potentials,  important  theoretical  and  experimental  methods 
using  the  J  integral  and  energy  release  rate  (originally  developed  for 
fracture  of  elastic  media  and  fracture  initiation  in  metals  with  plastic 
deformations)  may  be  extended  to  fracture  initiation  and  crack  propagation 
in  monolithic  and  composite  materials.  In  Section  2.2,  graduate  student 
Randy  Weatherby  describes  a  new  finite  element  model  for  analyzing  crack 
growth  in  materials  which  are  characterized  by  work  potentials.  It  is 
believed  that  both  the  use  of  the  "failure  zone"  in  a  finite  element  model 


and  the  study  of  path  independence  of  the  J  integral  with  crack  propagation 
are  new.  Graduate  student  Richard  Tonda  describes  in  Section  2.3  his  work 
on  determining  work  potentials  for  a  graphite/epoxy  composite  using 
unidirectional  bars  and  circular  tubes.  The  computer,  computer  programs, 
and  reduced  data  were  all  lost  in  a  fire  on  December  31,  1984.  This  study 
has  been  discontinued  for  now;  but  it  probably  will  be  restarted  early  this 
coming  Fall.  In  place  of  this  work,  we  have  begun  a  study  of  flat  angle- 
ply  laminated  bars  under  combined  axial  and  torsional  loading,  as  described 
in  Section  2.4.  A  third  graduate  student,  Mark  Lamborn,  recently  joined 
the  project,  and  he  has  contributed  to  the  study  described  in  Section  2.4. 
Preliminary  data  confirm  the  existence  of  a  work  potential.  Also,  early 
results  indicate  that  the  tests  of  flat  bars  with  axial  and  torsional 
loading  will  be  useful  for  Mode  III  delamination  investigations. 
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2.2  Analysis  of  Crack  Growth  in  Damaged  Media 

Using  a  Generalized  J-Integral 

The  purpose  of  this  research  is  to  develop  and  study  analytical 
methods  for  use  in  nonlinear  fracture  problems  with  distributed  damage.  We 
feel  that  these  methods  can  be  applied  to  a  wide  variety  of  materials 
including  metals,  polymers,  and  fibrous  composites.  Our  current  approach 
is  to  deal  with  crack  advance  purely  from  a  mechanics  viewpoint,  given 
certain  characteristics  of  the  zone  of  failing  material  at  the  crack  tip. 

A  central  idea  in  the  theory  of  fracture  which  we  follow  is  that  the 
stresses  and  strains  in  the  continuum  are  bounded  rather  than  singular  at 
the  crack  tip.  This  condition  is  met  by  the  introduction  of  a  "failure 
zone"  at  the  end  of  an  advancing  crack.  The  failure  zone  is  a  thin  layer 
of  highly  strained  and  damaged  material  which  is  usually  surrounded  by 
inelastic  material  (see  Figure  1).  The  failure  zone  can  be  represented  in 
a  finite  element  model  as  a  nonlinear  foundation  which  extends  ahead  of  the 
crack.  This  foundation  is  divided  into  discrete  elements,  which  we  call 
failure  zone  elements,  and  a  relationship  between  the  tractions  acting  on 
the  crack  surface  and  the  opening  displacement  is  specified.  The  traction- 
displacement  function  is  defined  so  that  after  a  certain  amount  of 
separation  between  the  crack  faces  the  surface  tractions  vanish  and  the 
crack  increases  in  length.  Similar  techniques  have  been  used  previously  in 
modeling  fracture  of  concrete,  adhesive  joints,  and  composites  [1,2].  In 
these  cases,  however,  the  continuum  surrounding  the  failure  zone  is  assumed 
to  be  linearly  elastic,  whereas  we  are  presently  considering  problems  in 
which  the  surrounding  continuum  is  a  rate  independent,  inelastic  material. 


Figure  1.  Crack  Tip  Region 


Under  certain  conditions  (which  will  be  stated  shortly)  it  can  be 
shown  that  there  is  a  path  independent  J-integral  [3],  and  it  is  equal  to 
the  total  work  input  per  unit  area  to  a  ligament  at  the  trailing  edge  of 
the  failure  zone.  First  consider  a  line  integral  defined  by: 


3u. 

JB  *  i  (tn‘-  Ti  H7>  ds  <11 

where, 

B  =  the  contour  of  integration 

$  =  potential  function  analogous  to  strain  energy  density. 

=  traction  vector  acting  on  B. 
u^  =  displacement  vector. 

n^  =  x^-component  of  the  unit  vector  normal  to  B. 

As  an  example,  using  j^fto  denote  the  value  of  (1)  for  the  closed  contour  in 
Figure  2,  and  expressing  it  as  the  sum  of  contributions  from  each  segment 
of  the  total  path,  we  have,  using  clockwise  integration. 


JAB  +  JBC  +  JCD  +  JDE  +  JEF  +  JFG  +  +  J’ 


GH  JHI 


JIJ  +  JJA 


(2) 


For  an  opening  mode  crack,  symmetry  can  be  used  to  give: 


2  <JAB  +  JBC  +  JCD  +  JDE  +  JJA^ 

For  the  material  model  we  use,  J^will  be  zero  when  evaluated  around  any 
closed  contour,  B,  provided  that  inside  of  B,  Tmax(x,y)  =  Tmax^  for 
points  where  t>t  max  if  Tmax>To'  here,  t  (x,y)  is  the  octahedral  shear 
stress,  t 0  is  the  value  of  t  at  which  appreciable  plastic  deformation 
occurs,  and  fmax  is  the  maximum  value  of  t  seen  by  a  material  point  over 


the  entire  history  of  loading.  It  is  assumed  that  the  regions  which 


crack 

surface 


violate  thej£=0  condition  have  been  excluded  from  the  interior  of  B; 
Equation  (3)  becomes,  after  setting  0, 

2  (J^g  +  JgQ  +JQJ  +  JDE)  =  -2Jja  (4) 

Further,  we  note  that  n^  =  0  and  =  0  on  segments  AB  and  CD.  When 
tractions  in  the  failure  zone  are  given  as  a  function  of  displacements,  it 
can  be  shown  that  [4] , 


2Jja  =  -2T 


where  21  is  the  work  per  unit  area  done  to  a  ligament  at  the  left  end  of 
the  failure  zone.  These  observations  allow  (4)  to  be  rewritten  as: 


2  (JBC  +  JDE^  =  2r 


In  the  remainder  of  this  report,  when  we  give  a  value  for  the  J- 
integral,  we  are  referring  to  the  value  of  the  line  integral  in  Eq.  (1) , 
but  integrated  only  over  an  open  contour  with  end  points  at  the  left  end  of 
the  failure  zone  (in  the  case  of  a  right  moving  crack  tip).  For  example, 
in  Fig.  2  the  value  of  J  to  which  we  refer  is  2  (JBC  +  JDE). 

In  order  to  verify  by  an  example  that  the  necessary  conditions  for 
path  independence  of  J  are  satisfied  after  crack  advance,  we  analyzed  the 
small  scale  yielding  problem  in  an  isotropic  material  using  the  finite 
element  method.  This  problem  also  serves  to  demonstrate  the  usefulness  of 
the  failure  zone  element  in  determining  crar:'  growth.  By  definition,  in 
the  small  scale  yielding  problem,  the  plastic  zone  size  is  small  when 
compared  to  the  crack  length  and  the  distance  to  the  nearest  boundary.  The 
displacement  field  far  away  from  the  crack  tip  is  specified  to  be  of  the 
form  found  for  a  crack  in  a  linearly  elastic  body  ie., 


Uj,  -  KjVr  ft  (©)  (7) 

Figure  3  shows  the  dimensions  and  boundary  conditions  used  in  the 
finite  element  model,  and  Figs.  4  and  5  give  the  continuum  and  failure  zone 
constitutive  functions,  respectively.  Substructuring  is  used  to  condense 
out  the  degrees  of  freedom  in  the  linearly  elastic  region,  leaving  only  the 
nodes  associated  with  the  elements  appearing  in  Figure  6.  All  continuum 
elements  have  eight  nodes,  and  displacement  continuity  between  incompatible 
elements  (i.e.,  where  3  elements  share  a  common  boundary  as  shown  in  Fig. 
6)  is  enforced  approximately  by  a  penalty  method.  A  J2  deformation  theory 
of  plasticity  with  elastic  unloading  is  taken  as  the  material  model  for  the 
continuum.  Figs.  4  and  5  show  the  various  material  constants  which  are 
required;  these  include  parameters  describing  the  uniaxial  stress-strain 
curve  for  the  continuum  and  the  traction-displacement  relationship  for  the 
failure  zone.  For  simplicity,  it  is  assumed  that  the  traction-displacement 
curve  is  independent  of  the  amount  of  crack  advance  (i.e.,  the  work  of 
fracture,  2r ,  and  crack  tip  opening  displacement  are  held  constant). 
Loading  is  accomplished  by  increasing  the  stress  intensity  factor.  This  is 
not  done  directly;  rather,  a  modified  Riks'  algorithm  [5]  is  used  in  order 
to  get  past  limit  points  due  to  the  finite  element  approximation  which 
sometimes  appear.  As  the  stress  intensity  factor  increases,  a  plastic  zone 
forms  around  the  crack  tip,  and  the  failure  zone  elements  separate  to  give 
crack  advance.  Figure  6  shows  the  area  where  significant  plastic 
deformation  has  taken  place  after  the  crack  has  grown  some  distance.  The 
term  "active  plastic  zone"  refers  to  the  material  points  currently  at 
yield,  and  "plastic  wake"  refers  to  the  points  which  have  yielded  and  then 
unloaded  as  the  crack  tip  passed  by. 

A  olot  of  the  dimensionless  elastic  "energy  release  rate" as  a  function 


of  the  dimensionless  crack  advance  is  shown  in  Fig.  7.  This  curve  is  often 
called  a  resistance  curve.  The  energy  release  rate,$ ,  is  defined  here  by 
the  equation: 

S-  K,2  <8> 

The  crack  advance  mentioned  here  is  the  increase  in  length,  referred  to  the 
right  end  of  the  failure  zone;  this  end  is  where  the  displacement  6  in 
Fig.  5  is  reached.  Results  are  shown  for  different  meshes,  one  containing 
15  elements  in  the  active  plastic  zone  and  the  other  containing  7  elements 
in  the  active  plastic  zone.  It  should  be  noted  that  while  $  increases 
as  the  crack  advances,  the  work  required  to  completely  separate  a  ligament 
of  material  in  the  failure  zone  is  constant.  The  increase  in  $  is  due  to 
the  formation  of  the  plastic  wake  which  tends  to  reduce  deformations  ahead 
of  the  crack  below  what  would  be  present  if  the  material  were  truly 
nonlinearly  elastic. 

Figure  3  is  a  plot  of  TmaxA0  as  a  function  of  x  for  three  values  of 
the  y-coordinate.  This  plot  shows  a  region  near  the  initial  crack  tip 
(x=0)  where  tmax  is  above  yield  and  also  is  a  function  of  the  x-coordinate. 
Since  this  violates  the  requirements  for  path  independence  of  J,  this 
region  must  lie  outside  of  the  contour  path  in  order  to  get  J  values  that 
are  path  independent  and  equal  to  2F.  The  area  of  path  dependence  is 
shaded  in  Fig.  6. 

When  the  J-integral  is  evaluated  on  a  contour  that  surrounds  the  crack 
tip  and  lies  in  the  linearly  elastic  region,  the  relationship  [4] 

-  1-v2  K  2 

w  E  KI  (9) 


should  hold  before  and  during  crack  growth.  This  relationship  is  found  to 


13 


be  essentially  satisfied  for  contour  1  in  Fig.  6  which  surrounds  the 
plastic  zone  and  wake  (cf.  Table  1).  The  values  of  J  coining  from  contours 
2  and  3  listed  in  Table  1  are  not  equal  to$;  however,  they  both  have 
approximately  the  same  value  and  are  equal  to  the  area  under  the  traction- 
displacement  curve  given  for  the  failure  zone.  As  expected  from  theory, 
contour  4  also  gives  approximately  the  same  value  for  J  as  contours  2  and 
3.  Figure  9  shows  the  relationship  between  contour  4  and  contour  1.  These 
results  confirm  the  path  independence  of  J  as  discussed  in  [3]. 

Future  work:  We  plan  to  study  how  parameters  related  to  the  failure 
zone  (eg.  om,  6m,  and  m  in  Fig.  5)  affect  the  resistance  curve  (Fig.  7)  for 
small  scale  yielding.  Also  to  be  studied  is  a  split  beam  problem  with 
plastic  deformation  in  which  the  plastic  zone  is  large  compared  to  the  beam 
height.  When  generalized  to  an  orthotropic  media  with  damage,  the  results 
for  the  split  beam  should  help  us  to  interpret  data  from  delamination  tests 
on  composite  materials. 
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Table  1. 

J-Integral  Calculations 
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Eq.  (9):  J  = 
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KI  =  1.32  x  2r 
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Contour  1  = 


E  through  K 


Contour  4  =  A  through  0 


Figure  9.  Distinction  between  contours  1  and  4 


2.3  Studies  of  Unidirectional  Bar  and  Angle-Ply  Tube  Specimens 

The  initial,  preliminary  phase  involved  testing  of  standard  uniaxial, 
30°  off-axis,  tensile  coupons  of  Hercules  AS4/3502  material  in  a  dry,  room- 
temperature  state  under  load  histories  which  included  periods  of  constant 
rate  loading  and  unloading  as  well  as  constant  load,  as  illustrated  on  page 
5  of  publication  No.  3  in  the  Appendix.  The  load  history  was  designed  to 
(i)  obtain  the  type  of  information  needed  to  determine  mathematical 
representation  of  the  "strain-energy-like"  potentials  for  loading  and 
unloading  up  to  the  stresses  causing  laminate  fracture,  and  (ii)  determine 
the  extent  of  any  viscoelastic  behavior  which  may  exist  with  constant  and 
changing  damage,  and  also  characterize  it  mathematically  if  it  is  found  to 
be  significant. 

These  preliminary  tests  indicated  that  the  viscoelastic  effects  (apart 
from  damage  growth)  are  measurable.  The  relative  magnitude  of  these 
effects  were  not  influenced  by  ply-thickness  variations  from  6  to  12  plys, 
nor  by  the  off-axis  angle  of  the  specimen,  varying  from  lcf  to  90°.  The 
periods  of  constant  load  in  the  test  histories  used  allowed  for 
determination  of  the  significance  of  the  viscoelastic  effects  in  examining 
the  differences  between  loading  and  unloading  stress-strain  curves.  The 
following  observations  were  consistent  across  the  sample  spectrum  involving 
different  thickness  and  fiber  angles,  6  H0.° 

(i)  At  50%  of  ultimate  levels,  between  50%  and  80%  of  the  observed 
deviation  between  loading  and  unloading  stress-strain  curves  was 
attributable  to  normal  viscoelastic  effects.  The  percentage  of 
viscoelastic  effects  tended  to  increase  as  the  off-axis  angle 


increased. 


(li)  At  75%  of  ultimate  stress  level,  between  30%  and  50%  of  the 
observed  deviation  between  loading  and  unloading  stress-strain 
curves  was  attributable  to  the  viscoelastic  effects,  and  again 
this  effect  tended  to  increase  as  the  off-axis  angle  increased. 

The  amount  of  viscoelastic  behavior  relative  to  damage  effects  made  it 
necessary  to  develop  a  rather  elaborate  computer-based  data  reduction  and 
analysis  scheme  to  quantitatively  characterize  the  deformation  response. 
This  scheme  uses  so-called  "pseudo-variables"  (in  this  case,  axial  pseudo- 
strains1')  which  is  computed  using  a  convolution  integral  of  the  form  [3J, 


=  v1  £  >  a!  dt' 


where  E(t)  is  the  uniaxial  relaxation  modulus  and  ER  is  a  constant, 
arbitrarily  selected  reference  modulus.  With  this  approach,  the  material 
was  then  characterized  precisely  like  an  elastic  material  with  damage,  by 
examining  the  relationship  between  stress  (a)  and  pseudo-strain  (eR).  The 
technique  is  believed  to  essentially  "filter  out"  the  viscoelasticity  at 
fixed  damage  levels.  The  necessary  software  has  been  developed  and  was 
used  during  this  period  to  begin  analysis  of  the  data  from  off-axis 
specimens  using  fiber  angles  of  10,  15,  30,  45,  60,  and  90  degrees. 

Around  November  of  1984,  we  began  consideration  of  the  second  phase  of 
data  collection.  It  was  considered  necessary  to  conduct  an  experimental 
and  analytical  study  of  angle-ply  laminates  under  multiaxial  loading  in 
order  to  verify  the  existence  and  use  of  a  work  potential  4>  for  modelling 
realistically  large  amounts  of  distributed  damage.  This  type  of  laminate 
undergoes  significantly  more  damage  prior  to  global  fracture  than 
unidirectional  samples,  and  therefore  provides  data  that  enable  a  more 
critical  evaluation  of  the  work  potential  theory.  The  unidirectional  ply 


studies  provide  baseline  elastic  and  viscoelastic  property  data  in  the 
absence  of  the  significant  residual  stresses  which  exist  in  angle-ply 
laminates,  and  thus  can  be  used  to  predict  baseline  angle-ply  response  with 
no  significant  damage.  Hercules,  Inc.  of  Magna,  Utah,  has  provided  four 
(4)  tube  specimens  of  angle-ply  laminates.  These  were  to  be  tested  under 
axial  load  and  torsion  using  the  MTS  hydraulic  axial-torsion  machine  at 
Texas  A&M.  If  Eq.  (12)  in  [6]  was  found  to  be  true  using  data  from  the 
tubes,  then  a  potential  was  to  be  constructed  and  expressed  as  a  function 
of  laminate  thickness-averaged  strains  and  the  extent  of  damage  using 
suitably  defined  damage  parameters  [3].  Physical  inspection  for  damage 
using  x-rays  and  ultrasound  were  to  be  carried  out  when  possible,  in  order 
to  correlate  physical  damage  stress  with  characteristics  of  the  stress- 
strain  behavior. 

The  first  set  of  tubes  to  be  studied  consists  of  two  t±30]2s  and  two 
[±60]  2s  tubes,  all  2  inch  dia  x  17  inch  long.  They  were  to  be  tested  and 
the  findings  analyzed.  As  the  initial  task  in  this  study,  it  was  decided 
to  acquire  and  validate  a  finite-element  routine  which  would  allow  for 
careful  and  optimized  usage  of  the  four  tubes  available,  since  the  very 
high  cost  and  long  lead  times  typically  required  to  obtain  additional  tubes 
made  the  tubes  quite  valuable  structures.  Analysis  of  the  uniaxial  data 
proceeded,  but  was  secondary  to  preparing  for  tube  testing  which  was  to 
begin  in  January  of  1985. 

On  the  morning  of  December  31,  1984,  during  an  ice  storm,  the  buildina 
at  the  Texas  A&M  University  Research  and  Extension  Center  (which  housed  our 
HP  1000  computer  being  used  to  analyze  all  the  data  and  to  run  the  finite- 
element  model  previously  mentioned)  burned  to  the  ground  in  less  than  one 
hour.  For  all  practical  purposes,  all  of  the  data,  computer  code, 


intermediate  results  and  back  ups  (which  were  stored  in  another  room  of  the 
same  building)  were  reduced  to  ashes  and  were  totally  irretrievable. 
Needless  to  say,  this  tragedy  dealt  a  severe  blow  to  project  progress.  The 
code  developer  was  unique  to  the  particular  advantages  of  the  HP  1000 
vector  processing  system  as  was  the  finite-element  code.  Subsequent  to  the 
fire,  attempts  have  been  made  to  regroup,  and  some  tube  testing  of 
isotropic  (aluminum)  tubes  was  carried  out  to  verify  the  test 
configurations,  data  acquisition  and  data  analysis  techniques.  Tube 
testing  could  then  proceed  if  a  suitable  simulation  and  finite-element  code 
could  be  developed.  Since,  however,  it  would  take  considerable  time  to 
transfer  all  the  code  and  existing  data  and  techniques  to  another  machine 
and  since  Texas  A&M  University  intends  to  relace  the  HP  1000  no  later  than 
September  of  1985,  the  proposed  plan  is  to  place  this  portion  of  the 
project  on  hold  until  the  new  processor  is  available. 
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2.4  Studies  of  Angle-Ply  Bars  under  Axial  and  Torsional  Loading 


Preliminary  studies  of  narrow  rectangular  plates  under  axial  and 
torsional  loading  have  indicated  that  they  may  be  very  useful  for 
investigations  of  work  potentials  with  microcracking  and  delamination 
fracture  behavior.  Specimens  with  different  fiber  angles  and  dimensions 
are  under  study.  Figures  10  and  11  give  some  data  from  plates  with  nominal 
dimensions  of  2.5M  long  x  0.25"  wide  x  0.13"  thick;  the  "nominal"  stresses 
and  strains  shown  on  the  axes  are  defined  below.  Plates  which  are  wider 
and  longer  and  have  other  fiber  angles  are  also  under  study.  It  should  be 
emphasized  that  even  though  there  are  stress  gradients  and  consequent 
nonuniform  damage  distributions,  the  proposed  theory  based  on  work 
potentials  may  be  experimentally  verified  and  then  employed  to  obtain  basic 
fracture  properties  from  such  tests  without  having  to  make  a  detailed 
stress  analysis  of  the  specimens.  However,  some  stress  analysis  will  be 
done  to  aid  the  determination  of  deformation  properties,  optimize 
specimen  dimensions  and  select  fiber  angles. 

An  important  part  of  the  initial  study  is  to  demonstrate 
experimentally  that  work  potentials  exist.  A  theoretical  basis  for  such 
potentials  is  given  in  [6];  but  it  is  based  on  certain  idealizations,  and 
therefore  it  is  important  to  provide  also  a  good  experimental  basis  for  the 
existence  of  potentials.  Considering  here  only  the  problem  of  elastic 
laminates  with  damage  under  increasing  axial  force  F  and  torque  T,  a 
necessary  and  sufficient  condition  for  a  work  potential  $  =  <J>(u,e)  to  exist 
(where  u  is  the  axial  displacement  and  9  is  the  rotation  about  the  axis 
between  the  grips)  is 
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ress-strain  curves  with  proportional 


Theoretical  prediction 


Representative  axial  stress-strain  curves  with  proportional  axial  and  torsional  straining 


if  F  =  F(u,0)  and  T=  T(u,6)  have  continuous  first  partial  derivatives  in 
their  arguments  [7,  p.  170].  It  is  assumed  there  are  at  least  limited  sets 
of  histories,  u (t)  and  6(t),  (such  as  proportional  deformations  wherein  u/6 
is  independent  of  time)  for  which  F  and  T  are  essentially  independent  of 
history. 

Before  using  Eq.  (11)  with  experimental  data,  it  is  helpful  to  replace 
the  variables  by  measures  of  stress  and  strain.  This  normalization  process 
eliminates  first-order  effects  of  specimen-to-specimen  dimension 
variations.  Specifically,  we  use  "nominal"  stresses  and  strains  defined  as 

a  B  F/bc  ,  t  =  3T/bc^  (12a) 

e  S  u/L  ,  y  =  ce/L  (12b) 

where  b  =  width,  c  =  thickness,  and  L  =  length  (between  grips).  For  the 
special  case  of  long,  thin,  homogeneous  specimens  (L»b>>c)  a  and  e  are  the 
axial  stress  and  strain  respectively,  and  t  and  y  are  the  in-plane  shear 
stress  and  strain  respectively  at  the  surface;  this  is  shown  in  [8]  for 
linear  isotropic  materials,  and  it  can  be  shown  co  apply  also  to 
orthotropic  materials  whose  material  axes  are  parallel  to  the  specimen 
edges.  The  variables  in  Eq.  (12)  are  useful  for  normalizing  data,  whether 
or  not  the  stated  conditions  apply.  Equation  (11)  becomes 

3o/Dy  =  3(t/3)/"'£  (13) 


This  equation  has  been  used  to  analyze  the  data  in  Figs.  10  and  11  by  first 
writing 


where  tq  =  tq(y)  is  the  shear  stress  for  e=0;  also  f=ffe,y)  in  which  f(o,Y) 
=0.  Next,  integration  of  Eq.  (13)  with  respect  to  y  yields 

a  =  ~  S  f(e,Y')dY"  (15) 

3e  o  0 

where  a  =  aQ(e)  is  the  axial  stress  whenY=0.  Thus,  the  quantity 

3  Y 

A  a  =  ---  f  He,y')dy'  (16) 

oZ  0 

is  the  change  in  axial  stress  due  to  the  torque- induced  shear  strain. 

The  procedure  used  to  check  for  the  existence  of  a  work  potential  is 
to  cross-plot  the  data  in  Fig.  10  so  as  to  obtain  f  (which  is  the  change  in 
shear  stress  due  to  axial  strain)  as  a  function  of  Y,  for  fixed  values  of 
£,  and  then  predict  the  modification  to  axial  stress,  Eq.  (16). 
Considering  the  limited  amount  of  data  presently  available,  it  is  desirable 
to  curve  fit  the  data  to  analytical  expressions  to  aid  the  needed 
interpolations  and  extrapolations.  It  was  found  that 

,Y  '  k  2.55  (B+Cy) 

/  f (e,Y  )dY  =  Ay  £  (17) 

o 

where  A,B,  and  C  are  constants.  Using  this  expression  in  Eq.  (16)  yields 
the  change  in  axial  stress  due  to  torsion.  Only  for  e/y  =  0.92  is  there  a 
significant  effect  of  torsion  prior  to  fracture;  the  prediction  is  drawn  in 
Fig.  11.  The  agreement  between  theory  and  experiment  is  relatively  good. 

In  the  series  of  tests  shown  in  Figs.  10  and  11  there  is  only  one 
specimen  for  each  deformation  history,  and  thus  the  small  differences 
between  most  curves  in  Fig.  11  could  be  as  large  as  specimen- to-specimen 
differences.  Nevetheless,  it  is  encouraging  that  all  of  the  predictions 
from  Eq.  (16)  turn  out  to  be  of  the  same  order  as  the  observed  differences 
in  normal  stress.  Although  not  needed  to  check  for  the  existence  of  a 
potential,  it  is  of  interest  to  note  that  when  there  is  little  or  no 
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coupling  between  torque  and  axial  load,  the  stress-strain  curves  obey  power 
laws  over  a  wide  range  of  strains,  as  shown  in  Fig.  12. 

Future  Work:  It  is  planned  to  continue  the  axial-torsion  testing  of 
angle-ply  bar  specimens  to  study  both  deformation  history  effects  and 
determine  the  range  of  histories  for  which  work  potentials  exist.  Mode  III 
delamination  will  be  investigated  also  using  specimens  with  and  without 
embedded  edge  delaminations.  In  most  tests  conducted  so  far,  fracture 
occurs  away  from  the  grips,  and  therefore  it  is  believed  this  study  will 
lead  to  basic  results  on  the  effect  of  globally  distributed  microcracks  in 
a  multiaxial  stress  state  on  macrocrack  growth  (e.g.,  the  edge 
delamination).  The  J  integral  will  be  used  with  the  data  to  characterize 
the  mode  III  fracture  initiation  and  propagation  behavior. 
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Abstract 

Methods  of  quasi-static  deformation  and  fracture  analysis  are  developed  for  a  class  of  nonlinear  viscoelastic 
media  and  sample  applications  arc  given.  Selection  of  the  class  of  media  is  guided  by  actual  rheological  behavior 
of  monolithic  and  composite  materials  as  well  as  the  need  for  simplicity  to  be  able  to  understand  the  effect  of 
primary  material  and  continuum  parameters  on  crack  growth  behavior.  First,  pertinent  aspects  of  J  integral  and 
energy  release  rate  theory  for  nonlinear  elastic  media  are  discussed.  Nonlinear  viscoelastic  constitutive  equations 
are  then  given,  and  correspondence  principles  w  hich  establish  a  simple  relationship  between  mechanical  states  of 
elastic  and  viscoelastic  media  are  developed.  These  principles  provide  the  basis  for  the  subsequent  extension  of  J 
integral  theory  to  crack  growth  in  viscoelastic  materials.  Fmphasis  is  on  predicting  mechanical  work  available  at 
the  crack  tip  for  initiation  and  continuation  of  growth;  some  examples  show  how  viscoelastic  properties  and  the  J 
integral  affect  growth  behavior.  Included  is  the  problem  of  a  crack  in  a  thin  layer  having  different  viscoelastic 
properties  than  the  surrounding  continuum.  The  Appendix  gives  an  apparently  new  constitutive  theory  for  elastic 
and  viscoelastic  materials  with  changing  microstructure  (e  g  distributed  damage)  and  indicates  the  conditions 
under  which  the  fracture  theors  in  the  body  of  the  paper  is  applicable. 


1.  Introduction 

Methods  for  characterizing  and  predicting  crack  growth  in  materials  which  arc  elastic 
(except  for  the  small-scale  inelastic  zone  at  crack  tips)  are  well-established  theoretically, 
and  considerable  experimental  confirmation  exists  (1.2].  The  methods  for  linear  media 
now  commonly  use  criteria  for  initiation  and  continuation  of  crack  growth  which  are 
expressed  in  terms  of  stress  intensity  factors  or  energy  release  rate,  f  or  nonlinear  media, 
especially  rubber  (2|,  energy  release  rate  is  often  employed. 

Fracture  theory  for  materials  exhibiting  large  scale  inelastic  behavior  is  considerably 
more  limited.  The  J  integral  theory  [1.3]  has  been  successfully  applied  to  initiation  of 
crack  growth  in  time-independent  (clastoplastic)  isotropic,  homogeneous  media  undet 
small  strains.  An  analogous  parameter,  the  C*  integral  [4],  has  served  in  a  similar  manner 
to  define  crack  speed  in  nonlinear  viscous  bodies.  Analytical  methods  and  their  verifica¬ 
tion  for  crack  growth  in  viscoelastic  media  are  mainly  limited  to  linear  isotropic 
homogeneous  materials,  although  some  theoretical  results  exist  for  linear  orthotropic  and 
nonhomogeneous  materials  [e.g.,  5  9|;  stress  intensity  factor  is  the  primary  characterizing 
parameter  for  fracture  initiation  time  and  crack  speed. 

The  objective  of  much  of  the  work  on  elastic  and  inelastic  materials  has  been  to  identify 
a  basic  crack-grow  th  controlling  parameter,  such  as  stress  intensity  factor  or  the  J  integral, 
which  accounts  entirely  for  the  geometry  of  a  body  (including  crack  geometry)  and  the 
applied  loads.  Values  of  the  parameter  which  produce  crack  initiation  and  various  crack 


speeds  are  then  found  for  a  given  material,  usually  experimentally.  This  information  max 
be  used  to  predict  crack  grow  th  in  different  geometries  and  aid  in  the  selection  and  design 
of  fracture-resistant  materials  and  structures.  It  is,  of  course,  very  important  th.it  the 
characterizing  parameter  account  fully  for  the  effect  of  geometry  and  loading  conditions 
on  crack  growth  if  empirical  corrections  for  each  application  are  to  be  avoided. 

In  this  paper  we  show  that  parameters  analogous  to  the  J  integral  and  energv  release 
rate  may  be  used  for  quasi-static  crack  growth  in  a  class  of  nonlinear  viscoelastic  materials 
under  finite  strain.  In  Section  2.  results  on  the  J  integral  and  energy  release  rate  for 
three-dimensional  deformation  of  nonlinear  elastic  media  are  collected.  Included  is  a 
simple  (and  apparently  new)  derivation  of  the  relationship  between  the  J  integral  and 
energy  release  rate.  The  J  integral  formulation  is  guided  by  Chen  and  Shield's  work 
[10,11],  and  interpretations  are  given  which  are  used  for  subsequent  application  of  •' 
theory  to  viscoelastic  fracture. 

Viscoelastic  constitutive  equations  and  methods  of  quasi-static  deformation  analxsis 
using  elastic  solutions  (correspondence  principles)  are  discussed  in  Sections  3  and  4. 
respectively.  Correspondence  principles  are  given  for  a  broader  class  of  problems  than 
considered  in  the  fracture  analysis;  for  example,  they  represent  a  new  approach  to 
analyzing  crack  closing  and  healing  phenomena  and  ablation  effects. 

Sections  3  and  4  provide  the  basis  for  using  the  J  integral  and  energy  release  rate  in 
nonlinear  viscoelasticity  problems.  This  generalization,  along  with  results  in  Section  5  for 
mechanical  work  input  to  the  crack  tip,  is  applied  in  Section  6  to  relate  fracture  initiation 
time  and  crack  speed  to  the  /  integral  and  viscoelastic  properties  of  the  continuum  arid 
failing  material  at  the  crack  tip;  these  relationships  represent  extensions  of  the  author's 
other  work  [12]  based  on  a  two-dimensional  J  integral  and  small  strains.  Then,  as  another 
application  of  the  theory,  in  Section  6  we  also  predict  the  effect  on  crack  speed  of  the 
rheological  properties  of  a  zone  of  damaged  or  otherwise  special  material  surrnundmn  the 
crack  tip.  In  practice,  viscoelastic  behavior  of  this  zone  is  often  significantly  different  from 
that  of  the  far  field  as  a  result  of  high  local  stresses,  dissipative  heating,  or  the  particular 
physical  situation;  the  craze  zone  near  crack  tips  in  glassy  polymers  [13]  and  an  adhesive 
interlayer  arc  important  examples.  Surface  or  path-independence  of  the  J  integral  exists 
for  materials  with  certain  types  of  distributed  damage  as  show  n  in  the  Appendix;  we  use 
this  important  property  here  in  accounting  for  behavior  of  the  material  surrounding  the 
crack  tip. 

The  deformation  and  crack  growth  theory  in  this  paper  is  not  much  more  involved  than 
that  of  nonlinear  elasticity  or  special  cases  of  linear  viscoelasticity.  This  simplicity, 
compared  to  what  one  would  expect  for  nonlinear  viscoelasticity,  is  a  direct  result  of  the 
particular  constitutive  equations  and  mechanical  variables  selected  to  characterize  rheo¬ 
logical  behavior.  We  believe  the  theory  provides  a  practical  approach  to  the  development 
of  realistic  damage  and  global  fracture  models  for  nonlinear  elastic,  viscous,  and  viscoelas¬ 
tic  media,  as  illustrated  by  the  author  for  a  particulate  composite  material  and  polycrystal¬ 
line  metal  [14.15]. 

2.  The  reference  elastic  problem 

Certain  basic  equations  for  clastic  materials  under  large  strains  are  summarized  in  this 
section.  They  are  expressed  in  terms  of  stresses  a(*  (t.y  =*  1,  2,  3)  and  displacements  uR 
referred  to  an  orthogonal  set  of  Cartesian  coordinates  \t  which  define  the  location  of 
material  points  in  the  undeformed  state  of  the  body.  (Although  B0  is  called  the 
“undeformed  state",  it  could  be  any  fixed  reference  configuration,  such  as  that  existing  at 
one  time  during  the  actual  deformation  history,  without  necessitating  a  change  in  the  basic 
theory.)  The  instantaneous  Cartesian  coordinates  »■*  of  material  points  in  the  deformed 
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body  BR  are  referred  to  the  same  fixed  axes  as  used  for  Bn,  and  therefore y,R  =  u,R  +  xr  In 
some  cases  o,R  or  uR  will  be  equal  to  the  stresses  or  displacements  in  a  viscoelastic  body; 
but  in  general  they  are  different,  and  therefore  the  superscript  R  is  used  to  make  this 
distinction.  Section  4  is  concerned  with  the  correspondence  between  states  of  elastic  and 
viscoelastic  bodies.  For  now  we  shall  just  list  and  discuss  relevant  equations  for  elastic 
media. 

The  stresses  aR  are  taken  to  be  the  components  of  the  so-called  Piola  stress  field  [1<S]. 
The  components  of  the  Lagrangian  stress  tensor  7jR  (17]  are  given  by  the  transpose  of  nR; 
viz.,  TR  =  oR.  Although  oR  is  not  in  general  symmetric,  these  components  are  very 
convenient  for  our  purposes  because  the  equilibrium  equations, 

3oR/3.v,+  7^  =  0  (!) 

and  the  relation  between  surface  tractions  7"R  and  stresses. 


are  identical  in  form  to  those  in  the  linear  theory.  (Throughout  this  paper  the  summation 
convention  is  followed  wherein  repeated  indices  imply  summation  over  their  range  unless 
stated  otherwise.)  All  quantities  are  referred  to  B„.  in  which  .v,  are  the  independent 
variables.  Namely,  F] R  is  body  force  per  unit  undcformed  volume.  7jR  is  the  surface  force 
per  unit  undeformed  area,  and  n:  is  the  outer  unit  normal  of  an  area  element  defined  in 
the  undeformed  state;  force  quantities  are  defined  as  the  vectors  existing  at  the  current 
time  /,  but  referred  to  the  undeformed  geometry.  B„. 

A  potential  0  is  assumed  to  exist  with  the  property  that 

or  =  30/3(mr,)  (3) 

where,  by  definition. 

3ii“/a.v;.  (4) 

For  an  elastic  material  0  is  the  strain  energy  per  unit  undcformed  volume  [17],  If  we 
invoke  the  physical  requirement  that  0  is  unaffected  by  rigid  body  rotation  (and  recall 
that  we  are  using  the  coordinates  .v,  as  the  independent  spacial  variables)  the  dependence 
of  0  on  the  displacement  derivatives  can  enter  only  through  the  symmetric  Green's  strain 
tensor  [26], 

£r  =  (1/2)[«r;  +  r/R,+  «R,«r.;]  (-M 

However,  considering  the  association  between  the  present  elastic  problem  and  the  actual 
(viscoelastic)  problem  introduced  later,  we  shall  not  restrict  0  hv  this  usual  physical 
condition.  Instead,  unless  stated  otherwise,  we  suppose  only  that 

0  =  0(  uR,.  v:,  .x ,.  r )  (6) 

implying  possible  dependence  on  the  nine  displacement  derivatives,  spacial  variables  x:. 
.v,  (allowing  for  nonhomogeneity  with  respect  to  ,x:  and  x , )  and  time  t  (allowing  for 
"aging"  changes).  The  body  is  assumed  to  be  homogeneous  with  respect  to  v,  for  now  to 
achieve  surface-independence  of  an  integral  that  is  useful  in  fracture  analysis.  With  the 
same  objective  in  mind,  we  assume  a  body  force  potential  0,  -  0,  («R.  ,x:.  ,v„  t)  exists, 
in  which 

f;r=  -  ;i0l/aMR  n) 

Consider  a  generic  volume  V  throughout  which  0  and  0,  exist  with  properties  defined 
by  F.qns.  (3)  and  (7).  Denote  the  bounding  surface  of  this  generic  volume  by  S  and  let  T,R 
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be  the  surface  tractions  and  nt  the  components  of  the  outer  unit  normal.  Multiple  (1)  by 
—  du*/d. integrate  over  the  volume  V.  and  then  use  the  divergence  theorem  (18]  to 
convert  the  integral  to  an  integral  over  the  surface  S'.  The  result  is 


d=j[(<t>- -  7]R3r^/a.v1]dr  =  0 


where  ds  is  an  area  element  in  the  undeformed  body.  When  <Fr  =  0  and  the  bod\  is 
homogeneous,  the  integral  ^/becomes  identical  to  the  v,  component  of  the  integral  vector 
derived  by  Knowles  and  Sternberg  [16];  however,  we  have  followed  Chen  and  Shield  [10] 
and  not  limited  <l>  to  dependence  on  strains.  By  considering  only  the  .x,  component, 
material  nonhomogeneity  with  respect  to  .x,  and  .t,  may  be  taken  into  account  while 
retaining  the  propertv  that  d=  0.  Only  this  component  is  needed  in  the  fracture  analysis  to 
follow. 

The  J  integral  and  crack  tip  model 

Suppose  the  body  contains  one  or  more  cracks.  Figure  1  shows  as  an  idealization  the 
intersection  of  a  crack  tip  region  and  local  crack  faces  with  the  plane  of  the  page  The 
dashed  line  is  the  intersection  of  a  representative  surface  i’  with  the  same  plane.  In  order 
to  meet  the  conditions  which  lead  to  d=  0.  no  crack  can  exist  inside  or  on  .S'. 

The  region  designated  as  the  failure  zone  in  Fig.  1  is  where  material  separation  occurs; 
it  may  contain  a  high  density  of  microcracks  or  microvoids.  The  material  comprising  this 
zone  in  the  undcformed  state  is  of  length  a  (not  necessarily  small)  and  is  assumed  to  exist 
in  a  layer  which  is  thin  (in  the  ,v2  direction)  relative  to  n.  Outside  of  the  failure  zone  it  is 
assumed  there  exists  at  least  a  small  neighborhood  around  the  crack  tip  for  which  d-  0.  A 
useful  definition  of  the  crack  tip  P,  and  one  that  we  employ,  is  that  it  is  the  leading  edge 
of  the  material  for  which  the  conditions  used  in  deriving  d-  0  are  not  met. 

It  should  be  noted  that  tractions  may  exist  along  the  crack  faces  to  the  left  of  the  failure 
zone;  for  example,  these  may  be  due  to  a  pressurized  fluid,  interfacial  friction  and  contact 
pressure,  or  damaged  material  connecting  the  faces  of  the  intact  continuum  (as  in  a  craze 
zone  in  some  plastics).  In  these  cases,  especially  the  last  one,  location  of  the  left  end  of  the 
failure  zone  (points  1  and  2)  is  somewhat  arbitrary  and  its  selection  may  depend  on  the 
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(a)  Undcformed  body.  IVf«»rmed  hods 

Figure  1.  Cross-section  of  crack  in  neighborhood  of  the  tip  P.  I  he  region  of  intense  damage  an  i  m.iiert. 

separation  processes  is  designated  the  failure  /one,  \shose  length  is  *\  in  ihe  undeformed  bod\ 
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particular  application  of  interest.  However,  normally  one  would  choose  it  so  that  a 
extends  at  least  over  the  part  of  the  surface  xz  =  0  for  which  crack  tip  material  behavior  is 
too  complex  to  be  able  to  predict  the  detailed  traction  distribution. 

Next,  we  introduce  the  integral  J,: 

r«  r.  3AuR 

<*> 

where  r,R  is  the  normal  stress  and  rR  and  tr  are  the  shearing  stresses  in  the  x]  and  v, 
directions,  respectively,  along  the  interface  between  the  failure  zone  and  continuum:  these 
are  Piola  stresses,  as  defined  previously.  It  is  assumed  that  the  failure  zone  is  sufficiently 
thin  that  the  stresses  in  Eqn.  (9)  are  the  same  along  both  top  and  bottom  portions  of  the 
interface.  The  AuR  are  the  components  of  the  relative  displacement  vector  between 
initially  adjacent  interface  points  across  the  local  crack  plane  x2  =  0.  The  viscoelastic- 
normal  stresses,  t2  and  a22.  and  displacement.  Am,/ 2.  along  the  upper  surface  of  the 
continuum  are  indicated  in  Fig.  2.  By  definition,  for  the  elastic  or  viscoelastic  problem,  the 
“relative  displacement  vector"  is  the  displacement  vector  at  the  failure  zone-continuum 
interface  above  the  local  crack  plane  minus  that  below  this  plane. 

For  later  use  we  shall  suppose  that  the  crack  tip  or  edge  P  is  essentially  straight  and 
parallel  to  the  .v,  axis  over  at  least  a  short  distance  /,  from  the  plane  of  the  page  (a,  =  0). 
Over  this  same  distance  it  is  assumed  the  failure  zone  integral  J,  does  not  vary.  Since 
n ,  =  0  along  the  top  and  bottom  interfaces  it  is  readily  shown  that  the  contribution  to  </. 
(8).  from  the  portion  of  .S'  along  the  top  and  bottom  boundaries  of  the  failure  zone  over 
the  crack  edge  of  length  /,  is  equal  to  -■//,.  assuming  the  integral  over  the  small  curved 
surface  at  the  tip  P  can  be  neglected:  this  latter  assumption  is  reasonable  as  long  as  the 
undeformed  failure  zone  layer  is  thin  relative  to  a  and  we  impose  the  physical  requirement 
of  finite  stresses  (including  those  at  P).  Thus,  from  the  condition  /=  0. 

./,  =7,  (10) 


y,  *  (i/;«)/v  +  -  ‘  R 


and  .S’,  is  the  portion  of  .S'  not  included  in  the  integration  along  the  failure  zone  over  the 
length  /,. 

Notice  that  y,  is  a  surface-independent  integral  in  that  its  value  is  the  same  (i.e.  Jt ) 
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regardless  of  the  choice  of  5,  except  for  the  conditions  stated  above.  It  can  be  reduced  to 
Rice’s  path-independent  /  integral  [3]  by  omitting  the  body  force  and  assuming  two-.' 
mensional  deformations.  Specifically,  let  S  be  a  cylinder  (having  the  cross-section  in  I 
1)  with  generators  and  with  normals  to  the  end  areas  which  are  parallel  to  x,.  At  the  ends 
«,  =  0,  and  therefore  the  contribution  to  /v  from  the  cylinder  ends  vanishes  if  7]R9wR/9.v, 
=  0:  this  condition  exists  on  the  ends  when  7]R  =  0  (e.g.  plane  stress)  or  7jR  =  TR  - 
3kr/9.X|  =  0  (e.g.  plane  strain)  or  7\R  =  9i/R/3.x,  =  9t/R/3.v,  =  0  (eg.  antiplane  strain). 
With  the  integrand  in  (11)  further  assumed  to  be  independent  of  a  ,. 

(4>  +  <J>.  )d*,-  rR“dA  (12) 

9.v, 


where  the  integration  path  C,  starts  at  point  1  in  Fig.  1  and  proceeds  counterclockwise  to 
point  2.  Assuming  0 ^ =  0,  crack  faces  parallel  to  .v,  and  traction-free,  small  strains  and 
rotations,  and  further  that  <J>  is  a  function  of  t/R,  through  the  strains.  (12)  reduces  to  the 
original  form  of  Rice’s  J  integral. 

Suppose  S',  in  the  three-dimensional  version  of/..  (11).  is  chosen  so  as  to  not  include 
any  portion  of  the  failure  zone-continuum  interface  (outside  of  /,).  We  may  then  consider 
(10)  to  give  a  basic  relationship  between  the  mechanical  state  of  the  continuum  through  /. 
and  the  characteristics  of  the  failing  material  along  a  segment  of  the  crack  edge.  In  some 
cases  one  may  want  to  use  a  failure  zone  integral  in  which  /,  includes  the  entire  crack  edge 
(or  a  large  segment  of  it).  If  the  integral  (9)  is  not  constant  or  the  edge  is  not  straight  along 
the  length  of  interest,  one  would  return  to  (8)  to  derive  the  desired  form,  accounting  as 
necessary  for  curvature  of  the  edge. 

Finally,  it  should  be  noted  that  (10)  and  (11)  do  not  depend  on  crack  face'  being 
parallel  to  the  a,  -  a,  plane.  Rather,  this  condition  is  imposed  only  on  the  laser  of 
material  comprising  the  undeformed  failure  zone. 


Energy  release  rate 

Up  to  now  we  have  not  considered  crack  growth.  By  introducing  a  virtual  crack  extension, 
the  value  of  /v  can  be  related  to  a  global  change  in  energy.  This  relationship  ma>  be  useful 
for  the  experimental  or  theoretical  determination  of  /.  for  elastic  and  viscoelastic  materi¬ 
als,  as  an  alternative  to  evaluating  it  directly  from  the  integral.  (11).  The  desired  equation 
may  be  derived  by  first  multiplying  (1)  by  a  change  in  displacement  St/,R.  integrating  over 
the  volume  V  of  B0.  and  using  (2),  (3),  (4),  and  (7)  along  with  the  divergence  theorem. 
There  results,  finally,  the  familiar  equation  for  virtual  work. 


Both  S  and  V  have  been  assumed  constant  in  deriving  this  result:  e.g.,  there  is  no  explicit 
change  in  crack  tip  location  or  phenomena  such  as  material  removal  through  melting. 
Crack  growth  will  be  simulated  through  a  suitable  choice  of  St/,,  as  in  [14],  thus  permitting 
the  use  of  (13).  An  edge  segment  of  length  (in  the  x,  direction)  of  only  one  of  possibly 
many  cracks  in  the  body  is  to  be  advanced  an  amount  Sa,  as  illustrated  in  Fig  3:  this 
advancement  is  assumed  independent  of  a,.  The  interface  between  the  failure  zone 
material  and  continuum  of  undeformed  planar  dimensions  (a  +  Set)  by  /,  is  denoted  as  .V, 
and  considered  to  be  a  portion  of  S  in  (13).  As  before,  we  denote  by  A mr  the  current 
relative  displacement  components  between  originally  adjacent  material  points  across  the 
local  crack  plane  (which  arc  specifically  the  displacements  between  the  portions  of  S. 
above  and  below  the  crack  plane).  Self-similar  crack  growth  is  now  imposed:  At/R  -  AnK(  £  ) 
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Recall  that  (10)  is  dependent  on  the  assumption  that  the  material  is  physically  homoge¬ 
neous  in  the  v,  direction  for  the  portion  of  the  continuum  bounded  by  the  failure  /one 
and  the  surface  S,  used  in  (11).  In  contrast,  this  assumption  is  not  needed  to  derive  (17).  in 
that  <1>.  <!>,  .  and  <t> ,  in  (18)  may  depend  explicitly  on  all  coordinates  .v,  (as  well  as  on 
displacement  derivatives  or  displacements,  as  previously  indicated). 

Finally,  it  is  to  be  observed  that  when  surface  tractions  on  .S',  or  bodv  forces  are 
specified  functions  of  v,.  or  treated  as  such,  the  potentials  are 

<J>,  =  -  <j>,  =  -  trmr.  ci  i 

liquation  (18)  then  takes  the  familiar  form. 

/>,  rKuH)de-j  7;VRdv.  ,::i 

bor  elastic  materials,  the  quantity  P,  is  the  potential  energy.  Its  physical  significance  is 
somewhat  different  for  the  class  of  viscoelastic  materials  discussed  in  the  next  section. 

3.  Viscoelastic  constitutive  equations 

The  constitutive  equations  which  will  be  used  are  based  on  (3).  but  the  displacements  i/H 
and  stresses  <rR  are  not  necessarily  physical  quantities  in  the  viscoelastic  body.  Instead, 
they  are  related  to  the  physical  displacements  u  (  \t.  t)  and  stresses  o,((.vA.  /)  through 
hereditary  integrals 

Specifically,  considering  displacements  first,  and  assuming  they  vanish  for  t  <  0. 

,,  ,  n  i)«, 

",  -*  f-H  J[  l-'('  -  T-  <)-jT -dr  (23) 

where  m,-«,(.v(.  t)  is  the  physical  displacement  in  terms  of  the  time  variable  of 
integration,  t,  and,  as  before,  the  coordinates  v,  of  the  undeformed  body.  The  quantity 
E  -  E(t  -  t.  t)  is  called  a  relaxation  modulus,  which  imparts  hereditary  characteristics  to 
the  deformation  behavior.  The  coefficient  /:R  is  a  free  constant  which  will  be  termed  t1 
reference  modulus;  it  is  helpful  in  discussing  special  material  behavior  and  introducing 
dimensionless  variables.  In  order  to  allow  for  the  possibility  of  a  discontinuous  change  in 
m,  with  time  at  t  =  0,  the  lower  integration  limit  in  (23)  and  succeeding  hereditary  integrals 
should  be  interpreted  as  0  unless  indicated  otherwise.  The  inverse  of  (23)  is 

r,  ,  3«R 

u,  -  Er)  D(t-T,t)~dT  (24) 

■'ll  VT 

where  D  =  Pit  -  r.  t)  is  termed  a  creep  compliance.  It  is  readily  shown  that  I:  and  D 
satisfy 

r„.T)dr-//fr-r„)  (25) 

J,„  >'r 

where  r„ /•  0  and  //(r -/„)  is  the  Heaviside  step  function-  //(/  0  and  1  fort-  /, 

and  t  >  /„,  respectively. 

It  will  be  helpful  to  use  abbreviated  notation  for  the  hereditary  integrals  For  any 
function  of  time.  /, 

\Edf)^Ey'  flAt  r.t)^c\r. 

{  /hi/)  V  Ey  fp(t  r.  r) '-{dr 

J.  ”  . 
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Also,  in  view  of  (23)  and  (24), 

/=  {  Ed{  Z>d/  } }  =  {  Z3d{  Ed/  } }  (26c) 

liquations  (23)  and  (24)  become,  respectively, 

uH  -----  {  Edi/,  ) ,  u,  =  {/>duK}  (27) 

Similar  hereditary  integrals  are  assumed  to  relate  a R  and  o;).  but  modulus  and 
compliance  are  interchanged, 

<»*  =  {  ;VK,  }•  °.,=  {A'tdoo}  (28) 

where  the  subscript  1  is  used  to  indicate  that  the  relaxation  modulus  and  creep  compliance 
(as  well  as  another  reference  modulus  £R1 )  are  not  necessarily  the  same  quantities  as  in 
(27).  Inasmuch  as  o;R  and  i/,R  are  not  in  general  the  physical  variables,  we  shall  call  them 
pseudo  stresses  and  pseudo  displacements,  respectively.  Similarly,  the  adjective  pseudo  will 
be  used  when  referring  to  the  potentials  <h,  (3).  <J>, .  (7).  and  <1^.  (19).  in  the  context  of 
viscoelastic  analysis. 

The  hereditary  integrals  used  here  are  linear  functionals  with  relaxation  and  creep 
functions  which  are  independent  of  .v,.  This  behavior  provides  the  useful  property  that 
differentiation  with  respect  to  a,  and  hereditary  integration  may  be  interchanged:  e.g.. 

flMR/il.v(=  (  Ed(ilu,/il.t, ) } .  a«„/a.v,  =  [  />d(i)uR/9.v  )}  (29) 

The  choice  of  constitutive  equation  (3).  with  pseudo  and  physical  variables  related  in 
accordance  with  (27)  and  (28),  is  motivated  by  the  fact  that  this  constitutive  theory 
approximates  well  the  deformation  behavior  of  various  materials,  and  leads  to  relatively 
simple  equations  for  viscoelastic  deformation  and  fracture  analysis.  The  latter  point  will 
be  brought  out  in  this  paper.  The  validity  of  the  constitutive  theory  has  been  discussed  in 
112.14.15)  for  the  case  in  which  o,R  =  an.  Observe  that  it  reflects  the  commonly  reported 
behavior  in  which  stress-independent  relaxation  or  creep  functions  in  single  integrals  serve 
to  characterize  Hereditary  phenomena  exhibited  by  many  nonlinear  materials.  Also,  some 
important  special  cases  may  be  readily  recovered  through  an  appropriate  choice  of  the 
material  functions.  Specifically,  if  E  =£',  =  />  1  =•  13,  1  -  —  /:R].  (27)  and  (28)  reduce 

to  u,  =  u  and  o;R  -  arl:  consequently.  (3)  reduces  to  that  for  nonlinear  elasticity.  If 

'  =  Id  =  ^R1  «nd 

D  --  ( /.  Er  )  '(/-  t)  (30) 

(where  r,  is  a  time  constant)  we  find  uR  -  o,.  and 

mr  tJUt,/"<d  (31) 

which,  together  with  (3),  yields  linear  or  nonlinear  viscous  behavior.  The  general  creep 
compliance  I)  -  D(t  t.  i ).  together  with  n^  =  all  (corresponding  to  /;,  -  ERX)  and  a 
pseudo  strain  energv  density  which  is  proportional  to  that  for  a  linear  elastic,  isotropic 
material  in  terms  of  u:Rr  yields  the  standard  constitutive  theory  for  an  aging  linear 
viscoelastic  material  with  constant  Poisson’s  ratio:  nonaging  behavior  results  if  fi  -■ DO  - 
t).  The  generalization  provided  by  (28)  in  which  pseudo  and  actual  stresses  are  not  equal 
is  useful  for  crack  closing  an  healing  analysis,  as  discussed  later 

It  should  be  mentioned  that  the  notation  D  -  D{t.  t)  is  employed  in  (14)  instead  ot 
/)  =  />(  t  -  r.  t).  These  forms  are  equivalent,  but  the  latter  is  more  convenient  in  the  studv 
of  crack  growth.  The  pseudo  energy  may  also  depend  explicitly  on  time  to  account  for 
effects  of  aging  in  the  nonlinear  behavior.  "Aging"  is  not  limited  to  intrinsic  material 
changes,  but  mav  be  due  to  direct  physical  causes  such  as  transient  tempeiatures  and 
residual  stresses  [  14). 


4.  Correspondence  principles 


Correspondence  principles  in  linear  viscoelasticity  theory  usually  refer  to  elastic-viscoelas¬ 
tic  relationships  involving  Laplace  transformed  stresses  and  displacements.  Instead,  here 
we  shall  give  three  correspondence  principles  for  time-dependent,  quasi-static  solutions  to 
nonlinear  elastic  and  viscoelastic  boundary  value  problems:  they  enable  a  viscoelastic 
solution  to  be  easily  constructed  from  an  elastic  solution.  In  terms  of  Piola  stresses  and 
coordinates  x ,  of  B0  the  equilibrium  equations  are 


au./a.v,  +  r,  =  o. 

(32) 

The  stress-pseudo  displacement  derivative  equations. 

a„=  {£,d[a<D/a(uR,)]} 

(33) 

and  body  forces 

(34) 

in  which 

t/R(  =  {  £d(3u,/a.v/)}.  «k  =  {  £d«, } 

(35) 

lead  to  three  intergro-differentia!  field  equations  for  the  three  displacements  u,  when 
substituted  into  (32).  The  functions  <!>  =  4>(m,r/,  xk.  I)  and  4>F  -  <!>,  (u,R.  xk.  i)  are  consid¬ 
ered  to  be  known;  until  Section  6,  where  the  Jv  integral  is  used,  we  allow  for  explicit 
dependence  on  all  three  coordinates  xk. 

As  boundary  conditions  we  assume  the  traction  potential  <!>,  =  ,v4.  i  )  is  speci¬ 

fied  on  a  portion  .V,  of  the  surface;  viz., 

«»„«,=  7>  -  {  £,d(9<l),/9MR)}  on  Sr  136) 

Elsewhere,  displacements  U,  =  U,(xr  I)  are  given. 

ti,~  U,  on  Sv.  (37) 

The  total  surface  is  S  =  Sf  +  Sk,.  Although  not  treated  here,  generalization  of  the  analysis 
is  easily  made  for  mixed  conditions  in  which  different  traction  and  displacement  compo¬ 
nents  are  specified  over  the  same  part  of  the  surface. 

In  all  three  correspondence  principles  the  reference  configurations  Bu  of  the  elastic  and 
viscoelastic  bodies  are  specified  to  be  identical  (with  identical  cracks,  if  any).  The  first 
correspondence  principle  is  restricted  to  time-independent  surfaces; 

CT-l.  The  viscoelastic  solution  (i.e..  the  stresses  and  displacements  in  the  viscoelastic  body 
which  satisfy  (32)  (37))  is 

o,,=  {  £,daR}.  u,  -  (  DduK  (  (3S) 

where  ok  and  m(r  satisfy  equations  of  the  reference  elastic  problem.  (1).  (3).  (4).  and  (7), 
together  with  the  boundary  conditions. 

r,k  =  -  9<frT/a«R  on  5,.  (34a) 

=  U,R  s  {  L'dL'  }  on  S\  .  (34b) 

It  is  seen  from  (39b)  that  we  first  transform  the  given  displacements  L[  (if  any)  usii. 
the  hereditary  integral  and  obtain  those  needed  in  the  elasticity  problem.  The  governing 
equations  of  elasticity  for  the  variables  ok  and  uR  are  then  solved.  This  solution  is  used  in 
(38)  to  obtain  the  viscoelastic  solution.  That  (38)  is  correct  is  easily  established  by 
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substituting  it  into  (32)-(37).  If  the  body  forces  /• ]  or  surface  tractions  7]  are  specified  in 
the  viscoelasticity  problem,  then  one  transforms  them  to  obtain  FR=\DidFl)  and 
T:R  =  {  />,d /; ).  after  which  the  elasticity  problem  is  solved  using  the  potentials  in  (21). 

When  S,  and  .S',  vary  with  time  certain  difficulties  arise.  If  3«,/3/  *■  0  on  ST.  then 
(59a)  does  not  result  in  the  correct  condition,  (36).  This  problem  may  be  seen  by  observing 
that  for  the  solution  in  (38). 

{  M®.1}}",*  {£idno"/}  on  5r  (40) 

and  therefore  *  7].  Other  difficulties  are  due  to  the  effect  of  past  values  of  7]  and  l  ] 
on  current  values  of  TR  and  i'R.  For  example,  consider  d.ST/dt  >  0  and  <!>,  in  (21).  (By 
this  shorthand  notation  we  mean  at  least  a  portion  of  SL  becomes  in  time  a  surface  on 
which  7]  is  given.)  Then  the  traction  T,R  cannot  be  predicted  for  all  t  >  0  from  the  given 
boundary  conditions  on  the  part  of  S  where  the  change  from  a  displacement  to  traction 
condition  occurs. 

The  next  correspondence  principle  is  for  this  type  of  boundary  value  problem,  but  we 
assume  £j  =  Hence,  {  £jd/ }  =/  for  all  functions  /.  and  therefore  the  constitutive 
equation  and  body  forces  are.  respectively, 

o„  =  a<i>/a(M*,).  -a4>,../a«f.  (4i) 

Also. 

7]  s  -  3<I>  ,/ch/J*  on  S, ,  w,  —  Ut  on  Sv.  (42) 

CP-ll.  if  dA,/dr  >  0.  the  solution  of  viscoelastic  equations  (32),  (35).  (41).  and  (42)  is 

u,  =  {  7>d i/R  }  (43) 

where  oR  and  uR  satisfy  the  equations  of  the  reference  elastic  problem,  (1).  (3).  (4).  and 
(7),  together  with  the  boundary  conditions  in  (39)  in  which  TR  =  T,. 

Verification  of  (43)  is  readily  accomplished  as  before  by  substituting  this  solution  into 
the  governing  viscoelasticity  equations.  Inasmuch  as  the  elastic  and  viscoelastic  stresses 
are  the  same  throughout  V  and  on  S  at  all  times,  no  basic  difficulties  arise  in  verifying  the 
solution  if  3n,/<)r  *  0  on  S,  and  in  determining  TR  when  (21)  is  used  for  <t>T.  However, 
the  present  class  of  problems  obviously  allows  for  crack  growth,  and  certain  physical 
questions  of  material  continuity  and  interference  have  to  be  addressed.  Pursuing  this 
point,  we  observe  that  the  relative  displacement  between  crack  faces.  Aw,,  in  the  viscoelas¬ 
tic  body  is  the  difference  of  displacements  in  (43)  evaluated  on  adjacent  crack  faces. 

Au,  =  {DdAHR}  (44) 

where  Au,R  is  the  displacement  difference  in  the  reference  clastic  problem.  Since  we  have 
specified  the  instantaneous  geometry  of  all  cracks  in  the  elastic  problem  to  be  the  same  as 
in  the  actual  viscoelastic  body.  Ah,  is  correctly  predicted  to  vanish  until  the  time  r, .  say. 
when  a  crack  tip  reaches  any  particular  physical  location;  this  follows  from  the  fact  that 
Ahr  --  0  at  this  same  location  when  t  <  /,  (assuming  prior  cracking  and  rejoining  of  the 
crack  faces  has  not  occurred)  which  in  turn  implies  the  hereditary  integral  in  (44)  vanishes 
when  /  <  /,. 

The  present  solution.  (43).  does  not  account  for  contact  or  rejoining  of  crack  faces. 
Rather,  it  could  predict  that  adjacent  crack  faces  pass  through  one  another  if  in  the  actual 
situation  they  rejoin  and  interfacial  compression  exists.  Following  such  rejoining,  the 
stresses  in  (43)  are  not  valid,  and  the  solution  may  become  much  more  involved.  We  shall 
not  consider  the  general  problem  here;  rather,  only  the  case  in  which  cracks  are  initially 
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open  and  then  close  or  shorten  through  a  healing  process  is  discussed.  I  he  third 
correspondence  principle  is  concerned  with  the  problem  d S ,  /d  r  s;  0.  and  it  is  limited  to 
the  case  E  -  ER :  < 32)— (37)  still  apply  except  i/f  =  ui  and  «R(  =  u  *  are  now  used. 


CT-lll.  If  d.9,/dr  <  0  and  if  r)n:/(h  =  0  on  .9, .  the  viscoelastic  solution  for  the  case  /.  =  EH 
is 


=  {  C,doR  }.  ",  =  "!* 


(45) 


where  oR  and  i/,R  satisfy  the  equations  of  the  reference  elastic  problem.  (1).  (3).  ( 4 1.  and 
(7).  together  with  the  traction  boundary  conditions  (39a)  and 


on  S. 


(4h 


Verification  of  (45)  is  made  as  before.  Observe  that  elastic  and  viscoelastic  displace¬ 
ments  are  now  equal,  while  the  stresses  depend  on  the  relaxation  modulus  /  Also,  it  is  of 
interest  to  observe  that  if  tractions  Tt  are  specified  on  crack  surfaces  (such  as  in  the  failure 
zone)  the  tractions  7]R  in  the  elastic  problem  are  different  from  the  actual  values  since 
rR={/),d7;}. 

The  Correspondence  Principles  II  and  III  are  not  limited  to  crack  problems,  l  or 
instance,  they  may  be  applied  to  problems  involv  ing  contact  between  different  continua  or 
ablation.  Furthermore,  even  though  they  are  based  on  apparently  different  constitutive 
equations  they  may  in  some  cases  be  used  for  the  same  material,  at  least  as  an 
approximation.  Indeed,  for  linear  viscoelastic  behavior  in  which  the  only  effect  of  aging  is 
in  the  relaxation  moduli  E  and  E,  (i.e.  <l>  does  not  depend  on  time  other  than  through  the 
displacement  derivatives),  the  two  constitutive  equations  are  easily  shown  to  be  equivalent 
if  £,  =  E. 

The  remainder  of  this  paper  is  concerned  with  crack  growth  analysis  for  materials 
obeying  (41 )  and  (42).  The  Correspondence  Principle  II  (CP-11)  and  J.  integral  theorv  w 
be  used  in  the  development  of  criteria  for  predicting  growth  initiation  time  and  crass 
speed. 


z.'. 


% 


5.  Work  input  to  the  crack  tip 

An  important  quantity  in  crack  growth  analysis  is  the  mechanical  work  available  from  the 
viscoelastic  continuum  for  producing  the  separation  (or  at  least  a  significant  change  of 
state)  of  material  in  the  failure  zone.  This  work.  IF’,,  will  be  defined  using  an  idealized 
model  of  the  crack  tip.  and  subsequently  expressed  in  terms  of  the  parameters and  Js  of 
the  reference  elastic  problem. 

Consider  as  before  a  slender  failure  zone  (in  which  a  is  large  compared  to  the  initial 
thickness  of  the  failure  zone  in  the  .v,  direction)  and  locally  two-dimensional  deformations 
plus  antiplane  shearing.  Before  discussing  IT,,  let  us  recall  that  J,  for  the  elastic  problem 
was  defined  through  a  line  integral  taken  along  the  instantaneous  interface  between  the 
continuum  and  failure  zone.  Furthermore,  the  surrounding  material  was  assumed  to  obey 
(3).  while  no  such  restriction  was  imposed  on  the  zone  itself.  For  the  present  purposes  of 
discussion,  let  us  suppose  this  zone  consists  of  the  thinnest  material  layer  lor  which  (3) 
provides  an  adequate  representation  of  the  surrounding  material.  Figure  4  depicts  the 
deformed  failure  zone  (in  the  opening  mode  for  simplicity)  using  a  solid  line  to  indicate 
the  interface. 

For  the  viscoelasticity  problem,  the  work  per  unit  undeformed  area  (in  the  v,  '. 
plane)  input  bv  the  continuum  to  a  given  material  element  (of  width  d  v,)  in  the  failure 
zone  from  the  time  the  crack  tip  arrives  at  the  element.  („.  to  the  time  the  left  end  ol  the 


j 
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figure  4.  Deformed  cross-section  of  crack  in  opening  mode  showing  current  interface  between  the  continuum 

and  all  elements  in  the  failure  /one  (  - ).  and  trajectory  of  the  interface  for  the  element  currently  at  the  left 

edge  of  this  /one  ( . ). 


failure  /.one  arrives.  is 


<’e 

"•■I 


()Aw, 


<1f 


df. 
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The  quantities  t  and  An,  (/  =  1.2.3)  are  Piola  stresses  and  relative  displacements, 
respectively,  along  the  interface,  and  therefore  they  have  the  same  significance  as  their 
counterparts  in  Jt,  (9).  Notice,  however,  that  differentiation  and  integration  in  (47)  is  for  a 
fixed  value  of  .v,,  whereas  that  in  (9)  is  for  a  fixed  time;  thus.  J,  does  not  in  general  reflect 
the  deformation  of  a  given  material  element.  As  a  result,  quite  apart  from  the  distinction 
between  elastic  and  viscoelastic  solutions.  IV,  is  believed  to  be  a  more  basic  parameter 
than  Jf  for  defining  material  failure. 

There  is  an  additional  important  difference  between  Wj  and  J,  as  they  have  been 
introduced.  Since  W'f  is  the  work  input  to  one  material  element,  the  continuum-failure 
zone  interface  has  to  consist  of  the  same  continuum  material  points  at  £  ■■■  0  as  at  £  -  a. 
The  dashed  line  in  Fig.  4  is  intended  to  represent  this  interface.  The  height  of  the  materia! 
element  at  £  =  0  is  indicated  by  /»„.  which  defines  the  thickness  of  the  layer  that  ultimatelv 
becomes  part  of  the  failure  zone;  depending  on  the  material  behavior,  the  two  interfaces 
may  essentially  coincide  along  a  partial  width  or  entire  width  0  <  £  ^  a.  or  may  not 
coincide  until  £  =  «  (/  =  t„).  Considering  the  analysis  to  follow,  in  which  Au,K  and  An,  are 
to  be  related  through  CP-11,  we  should  use  the  same  material  interface  in  the  elastic  and 
viscoelastic  problems;  the  outer  interface  (the  dashed  line  in  Fig.  4)  may  be  used  if  the 
undeformed  value  of  /;„  is  small  enough  to  not  invalidate  (10).  We  shall  assume  that  this  is 
indeed  the  case. 

According  to  CP-11,  the  relative  displacements  in  the  elastic  and  viscoelastic  problems 
satisfy  (44)  if  the  tractions  acting  on  the  elastic  and  viscoelastic  continua  are  the  same. 


(Recall  that  a  traction  potential  <t> ,  was  used  in  the  boundary  conditions  of  (  i’  ll  icf 
(42));  but  it  is  clearly  sufficient  to  use  the  same  tractions,  regardless  of  whether  or  not  '.lies 
are  actually  known  or  are  expressed  in  terms  of  a  potential.)  Noting  that  stresses  ~  .uni  r  * 
are  tractions  along  the  continuum  surface  (apart  from  a  sign  change  due  to  n , ).  we  speeds 
t,r  =  rt  in  order  to  be  able  to  use  (44)  in  the  next  section.  Criteria  for  the  time  /,  at  which 
crack  growth  initiates  and  for  the  speed  of  propagation  a  will  be  studied. 

6.  Analysis  of  crack  growl h 

Initiation  of  growth 

For  initiation  we  consider  a  material  element  which  is  at  the  left  end  of  the  failure  /one 
£  =  a  from  /  =  0  to  As  the  body  is  loaded  (beginning  at  t  =  0)  the  crack  tip  P  moves  to 
the  right  but  the  element  at  £  =  a  does  not  necessarily  break  immediately.  Rather,  it 
breaks  at  t  =  /,,  the  so-called  initiation  time.  We  are  interested  in  expressing  the  work 
input  to  this  end  element  in  terms  of  the  far-field  parameter  y, .  In  order  to  simplify  the 
analysis  so  that  viscoelastic  effects  may  be  shown  clearly,  stresses  in  the  material  in  the 
failure  zone  will  be  idealized  to  those  of  a  time-independent,  rigid-plastic  body;  viz.,  we 
assume  the  r,  are  independent  of  t  and  £.  Thus,  from  (9).  (44).  and  (47). 

J,  =  T,Au*a  (481 

and 

tr, -  t,Ah„  =  T,  {  ZMA«*  }  =  {  DdJ, }  (491 

where  AuR  and  A //,„  are  the  relative  displacements  at  £  =  «.  The  last  result  together  with 
Jy  =J„  (10),  and  the  notation  in  (26b),  yields  an  explicit  formula  for  work  input  at  t  =  /,; 

fFf=£Rjf'o(r,-T,r,)^dT.  ('<)) 

(Sec  [12]  for  more  general  models  of  failure  zone  behavior.) 

This  result  is  the  crack  tip  work  per  unit  undeformed  area  in  terms  of  two 
continuum-related  parameters,  /v  (which  accounts  for  the  geometry  of  //,,  and  applied 
loads)  and  £>,  the  creep  compliance  of  the  continuum.  The  y,  integral  is  the  same  as ./  for 
an  elastic  material  when  expressed  in  terms  of  the  surface  tractions.  Also,  just  as  for  an 
elastic  material,  if  the  failure  zone  size  «  is  small  compared  to  all  other  geometric  features. 
yv  is  essentially  independent  of  failure  zone  size  and  properties.  Recall  that  f'R  is  a  free 
constant,  and  may  be  selected  as  desired.  (Its  value  does  not  actually  affect  U’f  because ./, 
turns  out  to  be  inversely  proportional  to  £R.)  Thus,  if  the  continuum  is  elastic  with  a 
constant  compliance  D.  we  may  use  £R  =  D  1  and  obtain  the  familiar  result  IF,  =  J 
An  equation  for  predicting  t ,  is  obtained  by  introducing  the  work  2F,  required  to  fail  the 
element  at  £  =  «.  Thus  (50)  becomes 

2I>{Ddyv’}.  (51) 

The  factor  of  2  is  used  because  the  “fracture  initiation  energy",  I’,,  is  defined  like  a  surface 
energy,  counting  the  cross-sectional  area  of  each  side  of  the  failed  element  at  £  -  n  as  one 
unit  of  area.  This  energy  is  not  necessarily  a  constant,  even  when  the  stress  state  of  the 
crack  tip  is  the  same  for  all  conditions  of  interest;  as  the  failure  zone  displacement  history 
may  affect  F,.  Whether  1’  is  a  given  constant  or  depends  on  At/,,,  (51)  is  an  implicit 
equation  for  predicting  r,  in  terms  of  the  history  of./,.  On  the  other  hand,  one  could  use 
(51)  to  obtain  T,  from  tests  of  laboratory  specimens  by  detc  mining  the  value  of  the  light 
hand  side  of  (51)  when  the  crack  starts  to  grow  under  various  test  conditions,  including 
different  modes  of  crack  tip  deformation. 
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Finally,  we  observe  that  for  the  opening  mode  of  crack  tip  deformation  in  a  localh 
isotropic,  linear  viscoelastic  material  in  plane  strain. 

J,  =  (1  -,:)A',2/£r  (52) 

where  A',  is  the  stress  intensity  factor  and  v  is  the  Poisson's  ratio,  assumed  constant.  (This 
familiar  expression  for  linear  materials  may  be  derived  from  (9)  and  (10).  for  bounded 
stresses,  by  using  [5,  Part  I  (15)  and  (22)j  for  Kt  and  displacement.)  Substitution  of  (52) 
into  (51)  to  express  the  initiation  criterion  in  terms  of  A',  yields  the  author's  earlier  result 
for  linear  viscoelasticity  |5,  Part  II  (64)]. 

Crack  speed 


Predictions  of  IP]  for  the  next  case  with  d  >  0  is  facilitated  by  using  £  rather  than  t  as  the 
independent  variable.  Thus.  (47)  becomes 


0A«, 


d£ 


(5.1) 


where  At/,  =  £).  Liquation  (53)  will  be  expressed  in  terms  of  Jy  for  short-term 

steady  state  conditions.  Namely,  the  speed  d  and  failure  zone  length  «  are  assumed  to  be 
essentially  independent  of  time  during  a  generic  time  interval  n/d  for  which  the  crack  tip 
moves  a  distance  «;  during  this  same  interval  it  is  further  assumed  that  r,  and  Aw,  are 
essentially  independent  of  v,  (although  depending  on  £).  a  is  small  compared  to  the 
distance  to  other  geometric  features,  and  that  there  is  no  significant  change  in  D(t  -  t.  /) 
due  to  aging  (through  the  second  argument,  t). 

These  conditions,  together  with  the  observation  from  (44)  that  A  it,  and  A  m,r  arc  related 
in  the  same  way  as  for  a  linear  viscoelastic  nonaging  material,  lead  to  the  approximation 
[5.  Part  111. 

Am,=  A'K /)(/.( )Aw*  (54; 

where  t  ?  ki/d.  The  factor  A  is  a  very  weak  function  of  slope  n  9  log  D/d  log  /.  and  is 
practically  1/3  for  the  entire  range  of  slopes  (0  s:  n  <  1 )  encountered  in  practice.  Equation 
(54)  and  the  value  of  A  stem  from  the  smooth,  cusp-shaped  relative  displacement  Amr(£  ) 
predicted  for  a  linear  continuum  (with  hounded  crack  tip  stresses),  such  us  that  illustrated 
in  Fig.  2  for  the  opening  displacement.  Although  further  study  of  the  accuracy  of  (54) 
seems  warranted  for  nonlinear  continua.  it  is  likely  to  be  a  good  approximation  in  mans 
cases  in  view  of  the  insensitivity  of  Au  /Aw,(  to  the  detailed  behavior  of  A wR(£)  (5.  Part 
II);  indeed.  (54)  may  not  require  a  to  be  small  or  constant. 

Substituting  (54)  into  (53)  and  using  the  same  type  of  approximation  as  in  the  linear 
theory  [5.  Part  II).  which  does  not  require  t,  to  be  spacevvise  constant,  we  find 

r"  9Amr 

»;  =/V>(r.t,,)jf  r,--^d£  (55) 


where 

=  ka/d.  (56) 

The  integration  and  differentiation  in  (55)  is  for  v,  fixed,  while  that  in  (9)  is  for  /  fixed 
(i.e..  fixed  crack  tip  location);  however,  because  Aw1'  and  r,  are  independent  v,  (or  t)  for 
short-term  steady-state  growth,  these  integrals  are  equal.  Equation  (55)  thus  reduces  to 

i)Jt.  r'7' 

E  pon  equating  If]  to  the  wor'  required  for  failure  of  a  material  element.  21  .  and  using 


(10)  we  obtain  an  implicit  equation  for  a. 

21'  =  t)Jv.  (58) 

It  should  be  recalled  that  the  crack  speed  has  been  assumed  constant  for  only  the  generic 
period  a/d.  Consequently,  a  as  well  as  the  other  relevant  parameters  may  vary  over  much 
longer  periods  without  invalidating  (58).  Similar  to  the  initiation  problem,  the  linear 
theory  (5,  Part  II]  is  recovered  from  (58)  by  using  (52)  for  7, .  Also,  for  a  linear  viscoelastic- 
material  whose  creep  compliance  obeys  a  power  law  in  time,  the  factor  k  in  (56)  may  be 
expressed  exactly  in  terms  of  Gamma  functions  if  the  stress  in  the  failure  zone  is 
spacewise  constant  [27,  p.  118]. 

The  failure  zone  may  be  viscoelastic,  and  therefore  I'  could  depend  on  a  as  well  as  other 
local  parameters.  After  allowing  also  for  possible  dependence  of  a  on  T,  and  a.  for 
proportional  stressing  (58)  provides  implicitly  the  functional  relationship  d  =  ci(  Js  ).  as 
discussed  in  [12]  and  in  the  final  subsection.  The  effect  on  a  of  the  geometry  of  the 
undeformed  body  /?,,  (in  which  crack  lengths  vary  with  time)  and  applied  loads  is  entirely 
accounted  for  by  the  instantaneous  value  of  Jv.  In  principle,  a  detailed  model  of  the  failure 
/.one  could  provide  this  function;  but  if  the  effect  of  fundamental  material  parameters  is 
not  of  concern,  one  would  normally  determine  a(  J )  experimentally. 

Effect  of  a  process  cone  or  interlayer  on  crack  speed 

Figure  5  shows  a  failure  zone  within  a  “process  zone"  of  length  ft.  The  latter  zone  is 
introduced  in  part  to  account  explicitly  for  the  fact  that  with  some  materials,  especially 
plastics,  there  is  a  zone  around  the  crack  tip  which  has  different  viscoelastic  properties 
than  the  surrounding  continuum  [e.g..  13],  Also,  Fig.  5  can  be  interpreted  as  a  model  for 
crack  growth  in  an  adhesive  layer,  in  which  the  adjacent  continuum  represents  the  two 
elastic  or  viscoelastic  adherends.  Besides  these  uses,  we  may  consider  the  process  zone  to 
be  a  failure  zone  as  previously  defined  (but  now  with  some  limitations  on  its  constitutive 
properties),  and  thereby  obtain  detailed  information  on  f  from  the  subsequent  analysis. 

For  the  class  of  damageable  materials  introduced  in  the  Appendix,  a  pseudo  strain 
energy  density  exists  for  the  response  to  both  loading  and  unloading.  Additionally,  T,  is 
independent  of  path  under  many  conditions.  We  shall  assume  one  path-independent 


I  igure  5.  Cross- section  of  crack  in  undeformed  boctv  with  crack  Up  P  embedded  in  material  laser  ipr.KC'  /ore  l 
luting  deformation  char.iclerrliis  different  from  surrounding  continuum  the  failure  /one  is  nol  ncvcssu.lv 
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integral  exists  for  the  process  zone.  and  another  exists  for  the  surrounding  continuum.  The 
creep  compliances  and  't>  for  the  two  regions  may  he  different,  and  therefore  these 
integrals  are  not  necessarily  equal. 

Local  two-dimensional  deformations,  appropriate  to  the  use  of  (12)  for  the  process 
/one.  are  assumed.  Furthermore.  C,.  is  now  taken  to  be  the  dashed  line  in  Fig.  5.  which  is 
adjacent  to  the  continuum-process  zone  interface.  Subscripts  “a"  and  “b"  will  he  used 
when  necessary  to  distinguish  between  quantities  for  the  failure  zone  and  process  zone, 
respectively.  A  letter  subscript  will  not  be  used  with  parameters  of  the  continuum. 

The  process  zone  is  assumed  to  be  thin  in  the  x,  direction  (relative  to  ft )  so  that 
arguments  similar  to  those  given  previously  for  the  failure  zone  mav  be  emploved:  this 
slenderness  greatly  simplifies  the  analysis  while  the  essential  physical  features  are  retained. 
Thus.  (10)  is  used. 


J.=J,  i. 

where  .1.  is  in  (11)  and  (cf.  Fig.  5). 
rlt  itAi/f 


(59) 


(60) 


Similarly,  for  the./  integral  in  the  process  zone. 

J^-Ju  (01 ) 

w  here 


./,  - 


‘  H 


(62) 


It  is  important  to  recognize  that  in  using  (59)  (62)  we  are.  in  effect,  neglecting  the 
contribution  of  the  vertical  segments  at  P  and  Q.  Fig.  5.  to  the  integral  in  (8)  for  each 
contour.  The  other  simplification  employed  is  that  the  Piola  stresses  along  the  top 
horizontal  line  of  each  contour  equal  those  along  the  bottom.  The  same  simplifying 
features  are  assumed  to  apply  to ./. when  evaluated  on  C,:  thus.  (12)  yields 


/■ii  (IAmJI 

yW,T'i’“ird" w- 


(65) 


where  l„  is  the  integral  along  the  vertical  segment  below  point  1  and  above  point  2  in  Fig. 
5.  This  contribution  is  retained  even  though  the  integral  along  the  vertical  path  at  Q 
(within  the  process  zone)  is  neglected  T  he  difference  is  due  to  the  value  of  <!>,.  at  $  -  o. 
which  mav  be  very  large  as  a  result  of  damage  in  the  process  zone  (cf.  .Appendix),  while 
'!>„  at  i)  -  0  will  be  small  if  7  is  small:  the  contribution  of  the  integral  of 

r, n, */<»•' ,  on  the  vertical  lines  is  considered  to  be  negligible  due  to.  for  example, 
smallness  of  ,|  at  rj  -  0  and  |7],,|  at  £  -  a  or  the  slenderness  of  the  process  zone. 

Considering  the  fact  that  the  vertical  paths  at  Q  in  the  process  zone  and  the  continuum 
have  been  neglected  and  recalling  conditions  for  which  d  =  0.  we  are.  in  effect,  defining  Q 
to  be  close  to  or  ahead  of  the  leading  edge  of  major  softening  processes  and  to  straddle  t lie- 
location  where  the  creep  compliance  D  and  function  <t>  (for  the  continuum)  change  to  Ph 
and  <J>h.  respectively  (for  the  process  zone):  of  course,  if  I)  —  Du  and  =  <J>b.  the  latter 
constraint  is  not  involved.  The  distance  over  which  the  change  in  functions  occurs  should 
be  small  compared  to  /?  or  at  least  be  such  that  the  result  of  interest.  (75).  is  not  sensitive 
to  the  physical  location  selected  or  calculated  for  Q. 

Approximations  like  that  in  (54)  will  be  used  to  evaluate  the  viscoelastic  displacement. 
Thus,  for  the  continuum  just  outside  the  process  zone. 

±n,  IH  t.  / )  An",  t  A  ijz  d 


(64) 
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the  process  zone  just  inside  the  continuum. 
An,h  =  ■fcW  'h-  ')Amrh.  /„  =  A  hT)/a 


(65) 


and  the  process  zone  just  outside  the  failure  zone. 


=  £r £>h(6..  ')Amr.  r,  =  kj/a. 


(66) 


In  view  of  the  local  steady-state  assumption  implicit  in  these  approximations,  a  -  h:  also, 
because  £Rb  and  £R  are  free  constants,  we  have  used  £Rh  =  £R. 

Let  us  next  relate  to  yv.  This  may  be  done  by  using  the  continuity  of  displacements 
and  normal  and  shearing  stresses  across  the  continuum-process  zone  interface:  Am  =  Au,b 
and  r,  =  r(h.  The  first  condition  together  with  (64)  and  (65)  yields 


Am 


R 

/h 


D(l  0 


Amr. 


(67) 


Substitute  this  result  into  (63),  use  the  same  type  of  approximation  that  lead  from  (53)  to 
(55)  (but  now  with  the  D- ratio  in  (67)  entering  in  place  of  D  in  (54)),  and  then  employ  (59) 
and  (60).  There  results,  finally. 


■A 


i. 

Dhi'w  1 ) 


(68) 


where 

ip  =  A  (i/a .  (hB  =  A  hP/a  ( 69 ) 

and  approximately  A  =  Ab  =  1/3.  Additionally  (cf.  Appendix). 


(70) 


where  the  integral  is  taken  upward  along  the  vertical  lines  at  £  =  a;  the  basis  of  the 
physical  significance  of  4>h  we  have  /„  <  0,  where  /„  =  0  when  there  is  no  damage  in  the 
process  zone.  It  should  be  observed  that  we  have  neglected  the  contribution  from  the  bods 
force  potential  because  of  the  slenderness  of  the  process  zone. 

The  work  input  to  the  failure  zone  is  given  by  (57)  after  replacing  D  by  Db.  J,  by  Ju. 
(62),  and  then  using  (61)  and  (68), 


Wt  =  £„/>„(/'„.') 


o('V')  j 

^h(  V- ' ) 


(71) 


where 

i„  =  A.a/a.  A,  =  1/3.  (72) 

liquating  Hj  to  the  work  required  for  rupture  of  an  element  of  unit  area  in  the  failure 
zone.  21'.  we  find  the  implicit  equation  for  a. 

2Th=  £r£>(//(.  t)J,  (73> 


where 

F.RDh(i^.r).  (74) 

The  quantity  l'b  is  the  “fracture  energy"  of  the  process  zone;  it  is  introduced  in  order  to 
write  the  crack  speed  relation.  (73),  in  the  same  form  as  (58)  for  the  failure  zone  alone;  I',, 
consists  of  the  failure  zone  energy  !',  the  work  of  damage  (given  by  the  integral  in  (74)). 


21'h 


21' 


EM’n'O 


/«Mv, 
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and  the  “local"  creep  compliance  Dh  in  terms  of  two  local  times.  tn  and  th/l. 

As  a  result  of  damage  in  the  process  zone,  the  fracture  energy  T  may  be  negligible  (or. 
equivalently,  the  failure  zone  may  not  exist).  Equation  (74)  then  reduces  to 


•'a 


(75) 


which  is  to  be  substituted  into  (73)  to  obtain  the  equation  for  b(  =  d).  This  result 
corresponds  to  using  Jyh  =  0  in  (68). 

If  H'(  in  (71)  is  negative,  i.e.. 


f  < 


hd-v2 


(76) 


there  is  insufficient  mechanical  work  available  to  break  material  elements  in  the  process 
zone.  In  this  case,  the  crack  will  not  propagate  or.  at  least,  quasi-static  steady-state 
propagation  cannot  exist.  Even  if  Wf  >  0,  steady-state  propagation  may  not  exist  when 
(73)  predicts  d./v/da<0.  An  examination  of  (73)  indicates  this  situation  is  a  physical 
possibility  without  having  to  consider  dependence  of  a,  0,  I',  and  <!>,,  on  d.  For  example, 
suppose  /4>h  =  0,  E)£R  =  1  and  that  a.  /?,  and  T  are  independent  of  a.  Also,  let 
k  =  ka  -  kh  so  that  tfl  ~  ihp  -  k0/d  and  t„  =  kct/d.  Equations  (73)  and  (74)  yield 

21  77h(//).  t) 

'dJ^,JpS)~  "  (/7) 


It  is  helpful  to  express  the  derivative  of  Jy  using  logarithms. 
-  nh  +  =  d  log  Jjd  log  a 

where 

3  log  Dh(s,  t ) 

. 0  log  s  at  '  = 

and 


0  log  f\(s,  i ) 
a  log  s 


at  .v  =  alp/P 


(78) 


(79) 


(801 


are  logarithmic  slopes  of  the  creep  compliance  curve  (for  a  fixed  age  i)  at  a  creep  time  of 
s  -  tp  and  the  shorter  creep  time  v  —  at.,/0.  For  real  materials  (na.  nh)  >  0.  Also  n  ,  <  if 
s  is  sufficiently  small;  in  this  case,  d./v/d«  <  0.  For  many  materials  a’  long  creep 

times  v.  and  thus  dJ^/dd  >  0;  the  function  ./,(«)  is  therefore  predicted  to  have  a 
maximum  at  an  intermediate  crack  speed. 


Prediction  of  process  and  failure  zone  sizes 

The  lengths  p  and  a  which  appear  in  the  argument  of  creep  compliance  m  the  equations 
for  crack  speed  are  not  necessarily  constant.  In  fact  they  are  related  through  (59)  and  (61 ) 
respectively,  to  the  Jy  integrals  and  to  a  measure  of  the  stress  in  the  two  zones.  In  |12]  it 
was  noted  that  for  the  failure  zone  alone  the  condition  Jy  =  J,  leads  to  such  a  relationship, 
which  in  turn  reduces  to  that  derived  earlier  for  a  small-scale  failure  zone  in  linear 
viscoelastic  media  using  Barenblatt’s  condition  for  finite  crack  tip  stresses  [A  Part  1| 
Equations  (59)  and  (61)  are  for  slender  zones,  bui  0  and  «  may  be  large  and  the  material 
may  be  nonlinear,  viscoelastic,  and  anisotropic. 

Mechanical  state  solutions,  including  zone  size,  were  derived  in  [12]  for  a  failure  zone  ,n 
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a  power  law  material.  Here,  we  use  the  same  approach  to  derive  /<  and  n.  Consider,  for 
example,  the  use  of  (59)  and  the  power  law  potential  (a  so-called  homogeneous  function  of 
degree  M  +  1 ). 

)  He!” '  '<i>(  , )  (si  > 


to  obtain  an  expression  for  ft  The  result  is  (cf.  derivation  of  (39)  in  |12]). 

i  ,  i  w 


ft  = 


I 


./ 

A 


(S2) 


where  <rh  and  o,,  have  dimensions  of  stress  and  are  independent  of  v,.  but  mav  varv  with  i. 
a.  etc.:  thev  are  introduced  to  express  the  interface  stresses  -  and  potential  «J>  in  terms  of 
dimensionless  functions/,  and 

T,  =  »!./,-  <h  =  n„<h,/.  (S3) 

Also.  I  [denotes  absolute  value  and 

/,=  f'(a/li/i);)/id;  ( S4) 

where  j  ^  w/fl*  and  are  dimensionless  interface  displacements. 

X,  =  sign(nh)|o„/nlJl  (Ss) 

As  in  (12).  it  can  he  shown  that  under  certain  conditions  /,  is  a  constant  (apart  from 
parameters  which  appear  in  <l>„.  such  as  A/  and  aging  time  t);  although  linear  strain- 
displacement  relations  were  used  previously,  the  form  of  the  mechanical  state  solutions 
and  the  conditions  are  unaffected  by  the  magnitude  of  the  strains  if  the  power  law.  (XI ).  is 
applicable.  These  conditions  are:  (i)  the  “shape  factor"  for  interface  tractions  is  given  as 
/,=/('.  g, )  where  r  =  r),//?;  (ii)  the  continuum  is  locally  homogeneous  (i.e.  <1>,,  is 
independent  of  v,.  other  than  through  <UiK/d.\.):  (iii)  the  process  /one  si/e  ft  is  small 
compared  to  the  distance  to  geometric  features  outside  of  the  /one:  (iv)  the  crack  faces  arc 
locally  traction  free,  apart  from  the  interface  tractions  r. .  Without  further  analysis  one 
cannot  quantify  “local"  and  say  how  small  ft  must  be.  However,  it  may  be  necessary  for 
there  to  be  a  neighborhood  of  the  crack  tip  on  the  order  of  10  l()0/i  in  which  conditions 
(ii)  (iv)  are  met.  The  process  zone  must  be  small  enough  that  the  remote  stress  field  is 
essentially  the  singular  solution  a u  -  /•  w  u '  h.  where  r  is  the  distance  (in  the  v,  -  v; 
plane)  to  the  crack  tip.  Condition  (iv)  may  be  relaxed  to  allow  for  spacewise  uniform 
tractions  on  a  linear  continuum  and  spacewise  uniform  normal  tractions  on  an  incom¬ 
pressible.  nonlinear  continuum:  the  extension  is  achieved  through  superposition  of  a 
uniform  stress  field.  Also,  if  /  ----- /( r,  g , .  a.  J.,  n)  in  condition  (i).  then  /,  =  /,(<)../,.  n  ). 

The  analvsis  in  [12|  is  for  the  opening  mode  of  crack  tip  deformation.  However,  the 
form  of  the  results  is  the  same  for  other  modes,  including  mixed-mode  deformation.  One 
mav  then  argue  hv  analogv  that  /,  is  constant  under  the  same  conditions  as  stated  in  1 1  2]  if 
proportional  stressing  exists;  i.e..  if  the  ratio  of  stresses  in  the  remote  singtilat  stress  field 
is  independent  of  lime.  , 

For  uniaval  stress  o^-strain  «[',  behavior.  (XI)  implies  ntf,  -  .oh  n,  lie:  il 

M  <•  I.  a„  mav  he  interpreted  as  a  yield  stress:  if  A/  1.  n„  is  a  modulus,  l-quation  (X2)  in 
turn  provides  the  relation  between  o„.  a  measure  of  the  intrinsic  strength  of  the  ptocess 
zone  nh.  and  a  measure  of  the  external  loading  ./.  If  the  process  zone  is  also  a  power  law 
material  (in  the  neighborhood  of  the  failure  zone)  with  nonlinear  exponent  '/.  a  similar 
analvsis  vields 
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"hero  the  parameters  are  analogous  to  those  in  <X2).  If.  however,  there  is  considerable 
nonuniformly  distributed  damage  surrounding  the  failure  /one.  it  is  not  likely  that  a 
power  law  nonlinearity  w  ith  a  single  exponent  M  will  be  a  good  representation. 

■/,  its  a  characterizing  parameter 

C  onsider  again  a  power  law  nonlinear  material  with  a  small  scale  process  zone,  in  w  hich 
the  damage  and  failure  behavior  of  material  elements  is  unaffected  by  stresses  and 
deformations  prior  to  the  arrival  of  the  tip  Q.  Fig.  5.  Further,  assume  Js  and  />(  -  a )  are 
essentially  constant  during  the  time  ft/d  required  for  the  process  zone  to  propagate  its 
length.  I'he  state  of  stress  in  the  remote  continuum  where  the  aforementioned  singular 
solution  applies  is  determined  solely  by  the  current  value  of./,  for  proportional  stressing, 
just  as  for  an  elastic  material.  This  dependence  implies  ./.  is  the  only  remove  field 
parameter  which  affects  A  it*:  the  process  zone  characteristics  of  course  affect  Au,k.  From 
(64)  we  see  that  the  same  conclusion  applies  to  the  viscoelastic  displacement,  An,,  except 
crack  speed  now  appears.  The  damage  and  failure  parameters  a.  ft.  F,  and  /<J>b  in  (73)  and 
(74)  may  not  be  constant,  but  any  variation  will  be  due  to  7.  and  d  (e.g.  (82)).  apart  from 
aging  or  environmental  parameters  such  as  the  external  temperature.  (Local  temperature 
changes  due  to  mechanical  deformation  are  determined  similarly  bv  ./.  and  d.)  Equation 
(73)  therefore  serves  to  define  the  function  d  =  d(Jv),  indicating  that  7V  is  the  ‘characteriz¬ 
ing  parameter"  for  crack  speed.  In  principle,  this  function  could  be  obtained  experimen¬ 
tally  through  measurements  of  speed.  However,  by  introducing  specific  models  for 
behavior  of  the  process  and  failure  zones,  one  could  use  (73)  to  relate  crack  speed  to 
material  parameters  as  well  as./..  Elementary  examples  are  given  in  [12]  for  a  failure  zone 
in  a  continuum  obeying  a  power  law  with  respect  to  both  time  (through  the  creep 
compliance)  and  strain:  the  relation  d  is  derived,  where  k’  is  a  simple  function  of 
both  exponents  which  depends  on  characteristics  of  the  failure  zone.  The  assumption  of  an 
“elastic-like"  failure  zone  for  opening  mode  propagation  (F  and  r,  independent  of  speed) 
was  shown  in  (5,  Part  III]  to  provide  a  function  d  =  d(  J. )  which  agreed  well  with 
experimental  data  on  a  crosslinked  rubber:  in  this  study  there  was  no  process  zone  and  the 
continuum  was  linearly  viscoelastic  (cf.  (52)).  See  also  (7.20.21]  for  linear  behavior. 

The  J \  integral  may  serve  as  the  characterizing  parameter  for  initiation  time  or  crack 
speed  when  some  of  the  previously  stated  conditions  (e.g.  small-scale  crack-tip  zone)  are 
not  met.  depending  on  characteristics  of  the  process  and  failure  zones.  An  example  was 
given  earlier  for  initiation  time.  (51).  in  which  a  was  not  restricted  in  size.  However. 
further  experimental  and  theoretical  studies  are  needed  to  establish  the  necessary  condi¬ 
tions. 

7.  Concluding  remarks 

Methods  of  quasi-static  deformation  and  fracture  analysis  have  been  developed  for 
nonlinear  viscoelastic  media.  The  correspondence  principles  which  provide  the  basis  lor 
the  analvsis  are  not  limited  to  crack  growth:  they  apply  to  crack  dosing  and  healing  as 
well  as  to  other  types  of  problems  involving  ablation  and  interfacial  contact  and 
separation.  However,  only  crack  growth  examples  arc  given. 

(  onstitutire  equations 

Nonlinear  effects  in  F.qn.  (41)  for  stresses  arc  characterized  by  a  potential  <!'  which  is 
analogous  to  strain  energy  density  I  his  pseudo  energy  is  expressed  in  terms  ol  the  history 
of  displacement  derivatives  through  Fqn.  (23).  rather  than  the  history  of  Green's  strains. 
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and  therefore  material  objectivity  may  not  be  satisfied  when  large  deformations  exist:  i.e.. 
depending  on  the  deformation  history  and  material  type.  <t>  could  be  affected  bv  rigid 
rotations.  Material  objectivity  can  always  be  satisfied  when  linear  strain-displacement 
equations  are  applicable.  In  this  case,  one  would  express  <t>  as  a  function  of  displacement 
derivative  history  through  the  strain  history  [12]. 

A  single  hereditary  integral,  (23).  is  used  in  (41)  to  account  for  viscoelastic  effects,  and 
therefore  some  details  of  the  complex  stress-deformation  behavior  of  many  materials  may 
not  be  followed.  However,  the  theory  does  contain  general  material-objective  representa¬ 
tions  of  the  important  cases  of  nonlinear  elastic  and  viscous  media  under  small  or  large 
deformations  and  the  common  type  of  linear  viscoelastic  material  which  is  characterized 
by  one  independent  relaxation  or  creep  function.  In  certain  problems  of  crack  grow  th  in 
linear  viscoelastic  orthotropic  materials,  the  several  independent  creep  functions  combine 
into  just  one  function  for  predicting  load-displacement  response  of  the  crack  plane  [8]. 
This  feature  enables  us  to  generalize  (58)  and  (73)  for  crack  tip  work  by  simply  replacing 
D  w  ith  this  group  of  creep  compliances. 

For  large  deformations  of  viscous  materials,  the  current  geometry  would  be  considered 
the  "  undeformed"  state  in  order  to  recover  the  classical  constitutive  equations  [17];  the 
basic  expression  for  relating  crack  tip  and  far-field  behavior,  (10).  is  not  invalidated  in  this 
case  if  the  opening  displacement  Am,  along  the  failure  zone.  Fig.  1,  or  process  zone.  Fig.  5. 
is  small  compared  to  the  length  a  or  (9.  respectively.  On  the  other  hand,  this  condition  of  a 
slender  crack-tip  zone  in  the  current  geometry  is  not  needed  for  an  elastic  material. 

As  discussed  in  other  work  [12.14],  the  nonlinear  viscoelastic  constitutive  equations 
used  here  in  the  fracture  theory  may  be  written  in  the  form  of  a  special  type  of  a  so-called 
modified  superposition  principle  employed  successfully  w  ith  polymers  and  metals  [e  g.  15, 
19],  We  have  introduced  effects  of  aging  and  microstructural  changes  (e.g.  damage)  in  the 
standard  modified  superposition  principle.  This  “aging”  is  not  limited  to  independent 
physical  or  chemical  processes,  and  may  be  used  to  account  for  differences  in  nonlinear 
behavior  at  short  and  long  times  and.  as  shown  in  the  Appendix,  to  account  for  damage 
characterized  by  Lebesgue  norms  of  deformation-related  parameters. 

In  view  of  these  extensions  of  the  standard  single-integral  representation  for  hereditary 
behavior,  and  the  important  limiting  cases  contained  in  the  theory  (including  specific 
microcracking  models  and  deformation  theory  of  plasticity  with  elastic  unloading),  it  is 
believed  the  constitutive  equations  arc  sufficiently  general  to  account  for  the  primary 
features,  if  not  all  details  of  actual  deformation  behavior,  of  a  wide  variety  of  materials. 
Nevertheless,  considerable  additional  study  is  needed  to  establish  the  range  of  validity  of 
the  equations  for  different  materials  and  conditions.  For  example,  it  would  he  interesting 
to  determine  for  rubber  the  accuracy  of  the  stress-deformation  relation  in  (41).  which  is 
similar  to  the  theory  developed  and  successfully  applied  to  rubber  by  Tschoegl  and 
coworkers  [22], 

The  "pseudo  displacements",  u*.  which  appear  in  (41)  and  (43),  are  related  to  the 
physical  displacements  u,  through  the  hereditary  integral,  (23).  If  the  constitutive  theory  is 
valid  for  global  response,  the  material  behaves  overall  as  an  elastic  body  when  loads  or 
stresses  are  expressed  in  terms  of  pseudo  displacements;  as  illustrated  in  [15],  this  type  of 
behavior  for  nonaging  materials  with  constant  damage  may  be  easily  checked  by  convert¬ 
ing  experimentally  measured  displacements  to  pseudo  displacements  through  (23)  and 
then  examining  measured  load-pseudo  displacement  diagrams. 

Crack  growth 

Considering  the  complex  states  of  deformation  and  damage  around  real  crack  tip-, 
indirect  determination  of  specific  constitutive  equations  using  specimens  with  stationary 
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and  propagating  cracks  may  be  an  important  complement  to  studies  of  specimens  under 
homogeneous  deformation  histories.  This  determination  would  be  aided  by  the  simple 
relation  between  elastic  and  viscoelastic  displacements.  (43).  and  the  fact  that  (41)  for 
stresses,  including  the  generalization  in  the  Appendix,  leads  to  relatively  simple  equations 
for  crack  growth,  (51),  (58),  and  (73). 

The  quantity  Jv  in  these  equations  has  a  very  simple  meaning  in  certain  cases.  If 
hereditary  behavior  of  the  entire  continuum  outside  of  the  process  and  failure  zones  can 
be  represented  by  the  one  creep  compliance  D  (or  relaxation  modulus  £),  the  ./„  is  the 
same  as  the  familiar  J  integral  for  nonlinear  elastic  materials  and  C*  for  viscous  materials 
when  expressed  in  terms  of  externally  applied  loads  (rather  than  displacements).  The 
pseudo  potential  energy  /\,  (18),  is  related  similarly  to  the  potential  energy  for  elastic 
media.  If  «,  =  u,R.  except  possibly  for  a  small  neighborhood  or  crack  tips,  and  -  dP^/dA 
are  essentially  the  elastic  J  integral  and  energy  release  rate,  respectively.  The  ratio 
2Y/ERD  from  (58),  or  2rb/£RD,  (73),  then  appears  as  the  fracture  energy  for  an  elastic 
material  if  Jv  =  -  dPv/dA ;  this  “apparent  fracture  energy”  may  depend  on  crack  speed 
through  D,  T,  or  Tb.  For  the  rubber  studied  in  [5,  Part  III],  all  effects  of  speed  come 
through  the  nonaging  form  of  compliance,  D  =  D(f0).  For  a  viscous  body,  C*  is 
approximately  three  times  the  length-averaged  power  input  per  unit  area  to  the  failure 
zone  during  short-term  steady  state  propagation  if  we  use  rv  =  1  in  (30)  [12];  it  can  be 
shown  that  the  factor  is  exactly  three  for  a  linear  viscous  material  if  the  stress  in  the  failure 
zone  is  spacewise  constant. 

Direction  of  crack  growth 

Prediction  of  the  direction  of  crack  growth  has  not  yet  been  discussed  in  this  paper.  Many 
of  the  relationships  hold  whether  or  not  the  direction  changes,  but  the  problem  is  too 
complex  to  treat  here  in  any  detail.  We  would  only  suggest  a  possible  approach.  Referring 
to  Fig.  3.  suppose  for  purposes  of  discussion  the  local  coordinate  axis  is  fixed  and  various 
relative  orientations  0  for  the  continuum  are  considered.  The  actual  O' s  for  initiation  and 
continuation  of  growth  may  correspond  to  the  predicted  directions  for  which  t,  is  a 
minimum  and  d  is  a  maximum,  respectively.  These  criteria  automatically  account  for  local 
and  global  material  anisotropy  through  the  variation  of  values  of  the  material  and  loading 
parameters  (such  as  T.  a.  and  7V  in  (58))  with  respect  to  0.  Also,  these  proposed  criteria 
reduce  to  the  well-established  one  of  maximum  energy  release  rate  for  crack  growth  in  an 
elastic  isotropic  body.  The  equations  for  crack  speed  have  been  derived  under  the 
simplifying  condition  of  short-term  steady  state  behavior;  with  crack  tip  reorientation  it  is 
likely  that  the  equations  will  remain  valid  if  6  is  small  in  magnitude  and  essentially 
constant  for  an  amount  of  grow  th  equal  to  the  length  a  or  /?  of  the  crack  tip  zone. 

Crack  Up  models 

Emphasis  of  the  fracture  analysis  in  this  paper  has  been  on  predicting  the  mechanical 
work  available  at  the  crack  tip  for  initiation  and  continuation  of  growth.  The  right-hand 
side  of  (58)  and  (73)  is  this  work  at  the  failure  zone  (w'ithout  a  process  zone)  and  process 
zone  edges,  respectively.  By  assuming  the  theory  in  the  Appendix  is  valid  for  the  process 
zone,  we  have  obtained  some  information  on  how  the  creep  compliance  Db  and  pseudo 
energy  <J>b  of  the  process  zone  affect  the  required  work  21  b.  (74).  Viscoelastic  behavior  of 
the  embedded  failure  zone  in  Fig.  5  is  reflected  in  the  value  of  I  .  and  it  may  be  different 
from  that  of  the  process  zone  and  continuum.  Because  the  process  zone  is  slender,  we  were 
able  to  obtain  a  relatively  simple  relationship,  (68),  between  its  Jv  integral,  J^h.  and  for 
the  surrounding  continuum  (or  the  adherends.  in  the  case  the  process  zone  is  actually  a 
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thin  adhesive  interlayer).  Also,  it  should  he  observed  that  (71)  for  available  work  allows 
for  distributed  damage  outside  the  process  zone  as  long  as  the  theory  in  the  Appendix 
applies.  This  feature  is  important  for  many  materials,  especially  composites. 

With  the  dual  crack  tip  zones.  Fig.  5.  one  can  account  for  a  distinct,  complex  material 
separation  zone  (a)  within  a  relatively  well-defined  layer  of  damaged  material  (/l).  This 
geometry  may  be  a  realistic  model  of  the  delamination  tip  region  in  fibrous  composites 
with  rubber-toughened  matrices  {23.  Fig.  3]  and  a  cracked  craze  layer  in  plastics.  As  a 
special  case  the  failure  zone  could  be  omitted;  one  may  interpret  this  situation  as  the 
original  one  in  which  only  a  failure  zone  exists.  (58).  but  with  an  explicit  viscoelastic 
representation  for  the  fracture  energy.  As  shown  in  |5.  Part  III)  this  energy  for  rubber  may 
be  independent  of  crack  speed.  When  this  is  true,  the  molecular  theory  for  fracture  energx 
of  rubber  in  its  elastic  range  (2)  serves  to  relate  I',,  to  molecular  parameters,  and  (73) 
brings  in  the  only  viscoelastic  effects  through  /);  the  notch-tip  diameter  (  —  50  A )  used  in 
(2|  is  to  be  associated  with  the  process  zone  height  in  the  undeformed  state. 

It  is  believed  the  theory  in  this  paper  will  be  helpful  in  developing  detailed  crack  lip 
models  which  relate  growth  behavior  directly  to  local  physical,  chemical,  and  mechanical 
processes;  a  possible  general  approach  would  consist  of  using  the  available  crack  tip  work. 
(71).  in  a  local  nonequilibrium  thermodynamic  formulation  for  the  process  and  failure 
zones.  Much  of  the  published  work  on  crack  tip  models  employs  the  classical  singular 
solutions  for  the  local  mechanical  slate.  However,  whether  propagation  is  continuous  or 
occurs  in  steps  (e.g.  24,  25).  use  of  continuum  mechanics  and  thermodynamics  with 
bounded  stresses  should  lead  to  more  direct  relationships  between  basic  material  parame¬ 
ters  and  crack  growth,  as  illustrated  here. 

Acknowledgment 

The  a*.  is  grateful  to  the  U.S.  Air  Force  Office  of  Scientific  Research.  Office  ol 
Aerospace  Research  (with  Major  David  ("•••>«'>’ow  as  Proeram  Matv>"cr)  for  sponsoring  this 
research,  and  to  Mr.  J.R.  Weatherbv,  graduate  resc.i,k.n  assistant,  for  his  review  of  the 
manuscript  and  helpful  comments. 

References 

[1]  I).  Brock.  Elementary  Engineering  Fracture  Meehan  it  s.  Third  revised  cdn..  Martinus  Nijln-ff  Puhi-shcu 
(19X2). 

[2]  G.J.  Lake,  in  Progress  of  Rubber  Technology,  Applied  Science  Publishers  I. id.  (1983)  89-14  V 
(3|  J.R.  Rice,  Journal  of  Applied  Mechanics  35  (1968)  379-386. 

(4|  J.f).  Landes  and  J.A.  Beglov,  in  Mechanics  of  Crack  Grow  th,  AS’I  N1  S  1*1*  590.  American  Soviets  for  Jesting 
and  Materials  (1976)  128-148. 

(5J  R.A.  Sehapcry,  International  Journal  of  Fracture  11  (1975):  Part  I.  141-159,  Part  II.  369-3X8.  Pa  it  111 
549-562.  1 

[6 J  A. A.  Kaminskii.  Prtk ladnaya  Mckhamka  16  (1980)  (English  translation  (1981)  Plenum)  3-26. 

[7]  R.M.  Christensen,  Theory  of  Viscoelasticity,  An  Introduction ,  Second  edn..  Academic  Press  (19X2). 

[8]  G.S.  Brockwav  and  R.A.  Sehapcry.  Engineering  Fracture  Mechanics  10  (1978)  453-468. 

1 9J  R.A.  Sehapcry .  / nternat tonal  Journal  of  Frael are  14  (1978)  293-309. 

[10]  F.II.K.  Chen  and  R.T.  Shield.  Journal  of  Applied  Mathematics  and  Phwus  7  X  MP,  28  ( 1977)  i  72 
Jl  1]  R.T.  Shield,  in  Finite  Elasticity ,  R.S.  Rivlin  Pd..  The  American  Soctctv  of  Mechanical  P.nemecis  1 19  ",  1  -10. 
|12|  R.A.  Sehapcry.  in  Etn  vetopedm  of  Materials  Stience  and  Engineering .  Pcrgamon  Press  Ml)  <1IJX4)  Mo 
appear). 

1 1 3]  i;  Passaging  Polymer  23  (1982)  754*760 

(14]  R.A.  Sehapcry,  in  / f AS' /  Advances  in  Acrmpdtc  Structures  and  Material 5,  S.S  Wane  and  W  T  R •.■*)*  >n.  I  «!■»  . 

The  American  Society  of  Mechanical  Fngineeis  (1981)  6-2W 
1 15]  R.A,  Schaperv,  in  Proceedings  Ninth  I  S  National  C^ngrew  of  Applied  Mechanu  1  he  Aiueits  in  Soviets  «.[ 
Mechanical  l.ngineers  (198?)  237-245 


.  ■  «  t  .  -7.  •"«  i* .  -V»' * 


Correspondence  principles  and  a  J  integral  for  viscoelastic  media  219 

(16]  J.K.  Knowles  anil  H.  Sternberg.  An  hire  for  Rational  Mcchanu  v  and  Ana!\  ms  44  ( 1972)  1X7-21 1 

( 17]  y  l  ung,  b oundutions  of  Solid  Mechanics.  Prentice- 1  lull.  Inc.  (19pS). 

(IS]  Ml)  Greenberg.  Foundations  of  Applied  Mathematics.  Prentice- 1  fall.  Inc.  (197X). 

(19]  W.N.  Findley.  Journal  of  Applied  Mechanics  4X  (19X1 )  47-54 

(20]  \V.(i.  Knauss.  in  /V. for  mat  ion  and  Fracture  of  High  Polymer v.  H.  Henning  KauM.li.  John  A.  Hassell  and 
Robert  !.  Jaffee.  (ed.O.  Plenum  Press  (1974)  501-541. 

(21]  l  .N.  McC  artney.  International  Journal  of  Fracture  15  (1979)  31-40. 

(22|  N.W.  Tschoegl.  Polymer  20  ( 1979)  1365-1370. 

(23]  W.L.  Bradley  and  R.N.  C  ohen,  in  Pnnecdmgs  Fourth  International  Conference  on  Mechanical  Rchanor  of 
Materials.  3.  C'arlsson  and  N.(i.  Ohkon.  Pergamon  Press  (19K3)  >95  604. 

( 24 J  R.W.  Hert/berg  and  J.A.  Manson.  in  Fatigue  of  Engineering  Plastics.  Academic  Press  (19X0)  160. 

[25]  A.J.  Kin  loch.  Metal  Science  ( |9S0)  305*318. 

(26]  L.F..  Malvern.  Introduction  to  the  Mechanics  of  a  Continuous  Medium.  Prentice-hail.  Inc.  (1964). 

[27|  J.C».  Williams.  Fracture  Meehanus  of  Polymers.  Halstead  Press.  John  Wilev  &  Sons  (19X4) 

Appendix 

Fffci  t  of  time  -  carving  mu  rostruetural  changes 

Certain  important  inelastic  effects,  besides  those  represented  by  the  creep  compliance,  can  he  taken  into  account 
in  the  Jv  integral  theory.  These  additional  effects  are  associated  with  ‘'distributed  damage"  or.  what  ma\  be  a 
better  term.  *’ microstructural  changes".  Special  cases  are  microcracking,  dislocation  motion  and  generation,  hole 
growth,  and  breaking  of  entanglement  points  along  polymer  chains.  Here,  we  propose  a  theory  of  inelastic 
behavior,  allowing  for  large  deformations.  The  formulation  is  in  terms  of  the  reference  problem,  which  is  the 
reference  elastic  problem  with  damage;  the  viscoelastic  variables  are  still  related  to  those  for  the  reference 
problem  through  (27)  and  (2S).  For  lack  of  a  better  short  name,  we  are  using  "damage"  when  referring  to 
changes  m  the  microstructure  or  fabric  of  a  material.  However,  the  specific  “damaging  process’  does  not  have  to 
he  identified  here,  it  could  include  healing  (decrease  in  damage)  as  well  as  other  changes  which  are  beneficial  to 
structural  performance. 

(  onstitutire  theory 

It  is  assumed  on  phvsical  grounds  that  a  strain  energy  potential  exists  when  the  damage  is  constant  (for  instance, 
during  unloading  without  significant  interfaci.il  friction  following  development  of  microc racking  in  a  composite 
or  following  plastic  deformation  in  a  metal).  Denoting  this  potential  In  <I>‘  .  lor  reasons  to  be  given  we  assume  it 
has  the  form 

\ 

■t>'  •!•„(  ..  .  /  )  .  V  (Nil 

n  I 

The  enter  '!»„  through  the  functions  /,,  /„<  m,  , ».  as  a  special  » ase.  /„  max  be  a  strain  .  r  a  stiam  iiv  an.uit 

U'c  further  assume  that  all  of  the  fiiiutuuis  are  smooth  enough  to  permit  the  vurnus  mathematical  opeiatsms 
which  are  performed  I  he  «!'  s  mas  ilepend  on  x,  and  /.  as  indicated.  but  this  material  nonhomogeneitx  and  aging 
will  not  be  cxphulk  slumn  m  sitbsei|ucm  wotk  unless  needeil  for  cl. mix  Aha  throughout  the  Appendix 
summation  user  n  is  not  intended  unless  is  used 

Hv  deluiition.  all  quantifies  (which  will  be  i ailed  damage  p.namerer'\  are  c  mslain  dining  a  constant 
damage  process  In  tins  case  (3)  and  (A.l )  vicld  the  Pioln  stresses. 

,i.i.  . i<t>  ,w 

"  •  V  (Ml 

•*(».".)  •  , 

Bx  further  definition,  the  onlv  other  process  which  can  occur  is  called  a  damaging  protew  and  is  such  th  it 
/  /  tor  till  parameters  /tii;  at  the  end  of  tins  Appendix  a  simple  generalization  will  be  made  m  xxhteh 

I  *  /,,  for  some  of  the  damage  parameters.  It  is  assumed  that  the  stresses  are  continuous  in  the  displacement 
derivatives,  and  therefore  \  A  2)  defines  the  stresses  with  constant  or  varying  damage 

|  hi*  tot  mutation  tor  constant  damage  max  be  visualized  very  easily  by  consuleting  t  uniaxial  stress 
ft  i  rf*  >. strain  *  (  «»/#[*  .U,)  equation  written  tn  the  form  o  /U.<v  )  for  unloading  from  a  maximum  spam 
f  the  first-time  loading  curve  with  changing  damage  is  n  -  f(<.t).  where  «,  «  is  both  the  maximum  spam 

and  the  currcnl  strain  A  point  on  the  stress- strain  curve  for  loading  a  M f.r\  r-  uhs>  an  end  poini  on  an 
unloading  or  reloading  curve.  «  /t»  c .  i  at  *  * , 


V  V  S 
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A  potential  4>D  exists  during  a  damaging  process,  where 
oJ-aO^u",)  (A.3) 

because 

a</3("?.,)-a°£/»(<)-  (a  4) 

Equation  (A.4)  is  easily  verified  for  the  stresses  in  (A.2):  tn  evaluating  the  derivatives  of  stress,  one  must  first  set 
Fcn  =  F„  and  then  differentiate  both  arguments  of  each  4>„.  It  can  also  be  shown  that  if  at  least  one  of  the 
functions  <S>„  were  to  depend  on  more  than  one  damage  parameter,  (A.4)  would  not  necessarily  hold.  In  order  to 
fully  define  the  stresses  by  (3)  in  Section  2,  we  write  4>  -  4>c  when  the  damage  is  constant.  <t>  =  4>D  for  a 
damaging  process,  and  4>  =  d>c  -  iJ)D  at  the  transition  between  the  two  types  of  processes. 

As  an  aid  in  constructing  D  from  (A.2)  and  (A.3)  and  the  transition  continuity  condition,  and  in  providing  a 
physical  interpretation  of  this  damage  theory,  let  us  rewrite  <>c  in  a  different  but  equally  general  form. 
Specifically,  eliminate  <I>„  (n  =  1 . N)  from  (A.l)  in  favor  of  new  functions  P„  and  h„,  where 


=  PAK.rc„)+Jr'h„(F„)dF„ 

•'n 


p„(f.f*)  =  o.  (Ao) 

Notice  that  (A.6)  implies  P„(FC„,  Fcrt  )  =  0,  and  thus  the  integral  in  (A.5)  is  equal  to  d>„  at  the  transition  point. 
Fn  =  Fw  The  potential  for  the  constant  damage  process  may  be  written  in  the  form, 

H 

L  p.(p.pc)  <A7> 


*P  =  *o+Z,  rtm*AF.)dFm.  (AR) 

.-I1'® 

The  stresses  for  constant  damage  are 

(A. 9) 

"  a(^)  3(«r/) 

Requirements  of  stress  and  potential  continuity  at  the  transition  point,  F„  =  Fc„ ,  together  with  (A.3),  yield  for  the 
damaging  process, 


*n  =  *o  +  Z  /S„(F„')df; 


»(■!!,)  *  »(-r,) 


fc„(F.)-9P.(F..Fllt)/9F„  when  Ftn=F„.  (M2) 

Observe  that  <J>CD  »  <t>D  at  F„  -  F„.  Equations  (A  6)-(A.12)  constitute  the  damage  theory  which  will  be  used  in 
the  remainder  of  the  Appendix.  In  using  (3)  in  Section  2  with  the  damage  theory,  it  should  be  kept  in  mind  that 
4>  =  <t>r  for  constant  damage,  <t>  -  <J>°  for  a  damaging  process,  and  that  stresses  are  given  explicitly  hv  ( A  9)  and 
(A.1I). 

Let  us  next  consider  the  relationship  of  ♦D  and  <J>C  to  mechanical  work  and  dissipation.  Without  bod\  forces, 
the  mechanical  work  input  per  unit  initial  volume  lor  the  reference  problem  in  any  given  time  inters  .il  r,  s.  f  %  rh 
is 

^('..'h)”/'1«M".R,)/3']d'  (A1-) 

which  may  be  readily  established  by  means  of  the  divergence  theorem  Take  /,  -  0  and  assume  the  material  is 
initially  in  its  undeformed  state  («,R,  -  0)  and  that  it  docs  not  age  during  the  period  0  to  (i.e.  <}>  depends  on 
time  only  through  tiR( ).  Without  loss  in  physical  generality,  we  may  set  4>u  -  4>r>  ■*  0  at  /  -  0  (  onsider  next  an 
arbitrary  number  of  intervals  of  damaging  and  constant  damage  processes,  starting  with  a  damaging  process 


\7vX-4'-. 


' "■  i  ,  ►  ,  ;-  .*« *  vJ JM 


'  rot. 
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Substituting  (3)  into  t  A.l  3)  and  using  continuity  in  <t  at  the  transition  point  between  each  process. 

M/,n(0-'b>  =  ‘J,D('b)  (A.  14) 

when  fb  falls  within  a  damaging  process,  and 

w/,?(0.<b)  =  4>c(/b)  (A. 15) 

when  rh  falls  within  a  constant  damage  process.  Thus,  for  nonaging  elastic  materials  with  damage,  4>n  and 
have  a  very  simple  physical  meaning.  Notice  that  if  o,R(rb)  =  0  for  a  period  of  time  during  a  constant  damage 
process,  (A. 9)  gives  implicitly  the  residual  values  of  u,Rr  During  this  period  one  can  think  of  as  the  mechanical 
work  “dissipated”  due  to  damage  (in  the  sense  that  it  is  the  net  work  input  for  one  or  more  cycles  of  loading  and 
unloading).  If  healing  occurs  during  “damaging  processes",  the  work  input  could  be  negative.  For  uniaxial 
stress-strain  behavior,  one  may  view  4>rVrb)  as  the  area  under  a  tensile  stress-strain  curve  for  loading,  and 
•b‘  (rh)  as  this  area  less  that  under  the  curve  for  unloading. 

Pursuing  further  the  physical  significance  of  the  damage  theory,  we  find  that  a  special  microcracking  case  in 
[14,  (154|)  may  be  obtained  from  the  present  theory  by  setting 

oR  =  <r,,.  'v  =  '■  T  ,  =  Fm .  *„-F,. 

=  s(Fn)(F,- Fm).  *((•)  =  0.  (A. 16) 

The  damage  parameter  Fm  is  the  maximum  value  of  F(,  considering  the  entire  deformation  history  up  to  the 
current  time  F.quation  (A  12)  yields  h,(  F, )  =  -  g(F,).  Also,  we  find 

*c  =■[!-*(  ^m)][  -  *wl*  jf’l>  -  *<  )ldf.  (A. n> 

and 

<H,’  =  jf''[l-g(F1')]dF1'.  (A  18) 

The  stresses  for  a  constant  damage  process  are 

<A.1«) 

For  a  damaging  process  set  Fm  =  F,  in  (A  19).  The  function  g(  F, )  reflects  the  softening  effect  of  microcracks:  for 
no  damage  g  =  0.  and  for  complete  damage  (uniform  failure)  g  =  1.  Observe  from  (A. 19)  that  F,  is  the  pseudo 
strain  energy  density  for  an  undamaged  material.  It  should  also  be  added  that  g  may  vanish  over  an  F,  range. 
0  4  Ff  r;  Ft).  say.  1  o  this  case,  the  “damaging  process"  would  not  actually  produce  damage  until  tl.  ■  energy  /  n  is 
exceeded.  According  to  (A. 17)  and  (A. 18),  4>c  -  <t>l)  ■=  F,  for  0  y  Fr<  FD.  as  expected. 

A  damage  parameter  which  is  more  genera!  than  Fm  appears  in  (14,15).  It  is  derived  from  viscoelastic  crack 
growth  theory,  and  may  be  written  in  the  form  of  a  so-called  Lebcsgue  norm. 


If  p  =  oc.  (A. 20)  reduces  to  l.p  =  Fm;  this  case  leads  to  the  previous  theory.  ( A.l 7) -( A.  19).  On  the  other  hand,  if 
P  is  not  infinite  but  is  at  least  moderately  large  ( p  >  4),  and  if  the  deformation  history  over  a  period  0  r  K  r,n  is 
one  in  which  Fm  -  F] .  then  (15). 

l.p  =  AFjl1  r  (A  21) 

where  A  is  essentially  constant.  If  (A. 21 )  applies  up  to  t  -  i,„,  and  F,  <  Fm  for  t  >  f„..  then 


at  least  for  a  limited  period  of  time  beyond  rm.  The  damage  model  for  which  (A. 21 )  and  (A. 22)  applv  mas  be 
used  to  generalize  ( A.17)-(A.19):  viz.,  replace  g(Fm)  with  X<  Fmt'm/r)  and  g(  F, )  with  g(F,t'  r ).  Hence, 
microcracking  for  which  p  is  finite  leads  to  an  "aging"  elastic  material  with  damage 


The  7,  integral 

Dv  dividing  the  continuum  into  regions  of  varying  damage  (l.c.  where  a  damaging  process  exists)  and  constant 
damage,  we  may  easily  extend  the  theory  in  the  bodv  of  the  paper  to  allow  for  damage  if  the  displacement  field  is 
sufficiently  smooth.  In  effect,  the  material  may  he  treated  as  nonhomogencous  whether  the  nonhomogcneitx  is 
intrinsic  or  is  due  to  damage.  As  justification  for  this  extension,  consider  first  the  integral  .  (8).  We  will  allow 
for  explicit  dependence  on  s , ,  xv  and  /  in  the  potentials  4,n  and  /’„  in  (A. 7)  mid  (A  8)  (besides  dependence  sin 
*,!;  but  dependence  on  x,.  cither  explicitly  or  in  F(„.  is  excluded. 
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At  any  given  time,  let  5  *  .Vk  enclose  a  region  of  the  KhF.  undergoing  a  umsl.int  damage  pr.*cess  Clearly 
</==0  if  is  used  for  <t>.  If.  for  example,  we  were  to  use  <!>'  fnun  the  special  mu,  roc  rucking  nvdel  i  \  |7|,  it 
would  he  necessary  for  the  maximum  energy  /„,  to  be  independent  of  »t.  although  Fut  could  vary  with 
(distance  from  the  local  crack  plane)  and  t ,.  Now.  let  S  e  .Vj,  enclose  an  adjacent  region  undergoing  a  damaging 
process.  Here.  h[n  »  Fn%  and  thus*/  “  0  if  <t>n  is  used  for  4>.  even  though  ma>  vary  with  a,,  this  dependence  is 
through  Fn  =■  /^(m,r,  )  and  therefore  it  causes  no  problem  Next,  take  S  «  \  -*•  .V,,.  and  consider  the  interface 
where  .S'v  and  arc  adjoining  surfaces.  The  contribution  to</along  the  interface  from  hoth  \  and  Vfl  vanishes  if 
0.  and  4»f  are  continuous  across  the  interface;  this  follows  from  th-.  fact  that  n,  and  f  K  are 
continuous  in  magnitude  at  the  interface  hut  have  opposite  signs  on  .V  and  .V,,.  Consequently,  if  the  transition 
between  a  damaging  process  and  a  constant  damage  process  is  sufficiently  smooth  (X)  and  <  Hb  hold,  and 
surface-independence  of  J%.  (11).  and  path  independence  of  7V.  (12)  exist,  as  an  obvious  generalization.  .S'  max 
surround  an  arbitrary  number  of  connected  /ones  with  both  types  of  processes 

It  is  important  to  recall  that  *  0  is  required  in  each  constant  damage  princess  /one  For  example, 

suppose  the  body  is  loaded  and  a  damaging  process  /one  with  nonuniform  damage  is  produced  t  pon  partial 
unloading,  this  /one  may  become  a  constant  (with  respect  to  time)  damage  pr»xcss  /one  but  *F%„  >'»,  and 
therefore  </for  this  zone  would  not  in  general  vanish 

As  an  illustration  of  the  use  of  this  theory,  consider  the  crack  gross th  problem  of  Tig  *  treated  m  the  bods  of 
the  paper.  Locally  steady-state  crack  propagation  is  assumed.  We  suppose  that  a  given  material  element  in  the 
process  zone  b  undergoes  a  damaging  process  as  the  crack  tip  P  approaches  it.  vxhich  changes  to  a  constant 
damage  process  by  the  time  the  left  end  of  the  failure  /one.  £  ~  a,  arrives.  Nosv.  by  definition,  the  steads  state 
condition  is  one  in  which  each  clement  in  the  constant  damage  part  of  the  process  /one  has  djmage  parameters 
which  are  independent  of  v,;  e  g.,  for  the  simple  microcracking  nnHiel  dFM  /fl.r,  =  0.  Path  independence  in  the 
process  /.one  h  then  permits  us  to  evaluate  from  7vh  (cf  (61 ))  using  the  contour  C,  in  Fig  5.  The  result  for  7,h 
is  in  (68),  where  /rt  is  the  line  integral  of  -  <t>h  at  £  -  «.  (70).  The  potential  <t>h  is  actually  <t»J,  if  a  constant 
damage  process  exists  at  £  -  or  above  and  below  the  failure  /one.  In  view  of  (A  15).  the  integral.  -  is  the  net 
work  input  to  the  process  zone  per  unit  volume  integrated  from  the  bottom  to  the  top  of  the  pr»vess  /one  If 
there  is  no  damage  in  the  process  zone  (/*„-())  and  the  pseudo  strain  energy  density  vanishes  at  £  <*.  then 

l„  —  0.  However,  dissipation  would  still  exist  for  a  viscoelastic  material;  it  is  reflected  in  the  speed  dependent 
creep  compliances  which  appear  in  (71). 

Energy  release  rate 

F.quation  (17)  applies  with  damage,  even  if  the  homogeneity  condition  in  xx  needed  for  Js  is  not  satisfied  This 
generalization  may  be  shown  by  retracing  the  proof  without  damage,  but  using  <t><  and  <J>n  for  where 
appropriate.  Only  the  virtual  work  (13)  needs  to  be  examined  as  the  subsequent  steps.  (14)-(I6).  are  unaffected. 
Self-similar,  virtual  crack  advancement  is  imposed  through  application  of  appropriate  surface  tractions  t o  derive 
(17).  even  though  in  an  actual  crack  growth  process  this  type  of  advancement  may  not  occur. 

The  validity  of  (13)  may  be  established  using  the  same  procedure  as  for  (8)  with  damage.  Namely,  divide  the 
body  into  constant  damage  and  damaging  process  regions.  Inasmuch  as  a  change  in  the  displacement  field  occurs 
due  to  6a,  thin  layers  of  thickness  on  the  order  of  5 a  have  to  be  excluded;  these  are  the  lasers  in  which  the  type 
of  process  changes  during  the  crack  advancement.  However,  if  4>  and  5m,  arc  continuous  in  x,.  the  work  and 
pseudo  strain  energy  associated  with  these  layers  is  of  order  (5a):.  i.e.  thickness  v  5m,.  and  thus  the  lasers  dv 
not  affect  the  result  in  the  limit  5m  -*  0;  other  contributions  to  (13).  including  5IVR  in  (16).  arc  of  order  5 

Generalization  of  the  damage  model 

So  far  we  have  considered  only  the  case  in  which  the  potential  for  constant  damage  is  given  In  (A."i  and  all  of 
the  damage  parameters  satisfy  =  /•„  or  else  all  are  constant.  The  model  may  be  easily  generalized  bv  adding 
one  or  more  similar  groups  of  terms  to  <1><  ,  with  each  group  possibly  representing  a  different  phvsic.il 
mechanism.  The  regions  of  constant  and  varying  damage  for  one  group  need  not  be  at  the  sunc  h'calions  as 
those  for  another  group. 

A  constant  damage  potential  which  has  a  direct  effect  on  residual  stresses  due  to  damage  is 

,v:"> 

where  Or  -4>,(m,r  )  and  the  subscript  c  indicates  that  mr,  is  constant  during  a  constant  damage  poncss  I  he 
stresses  corresponding  to  this  potential  are  constant. 


(  d«t>, 

a, .  ~  I - 

\ 


'.'-V-v-v-v 
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For  a  damaging  process  u,R/c  =  «,R(.  and  a  potential  exists  whose  value  is  4> 11  ■=  <t>,(  uR, );  the  associated  stresses 
are  given  by  (A. 24).  but  without  the  subscript  c.  These  potentials  may  be  added  to  (A. 7)  and  (A. 10)  to  generalize 
the  damage  theory.  More  than  one  potential  4>r  may  be  used,  with  each  depending  on  different  combinations  of 
the  displacement  derivatives.  As  one  example,  we  obtain  the  Hencky  deformation  theory  [26J  with  linear  clastic 
unloading  if  (A. 23)  is  added  to  only  <t>w,  where  <t>„  is  taken  as  the  strain  energy  density  for  an  isotropic  linear 
elastic  material,  we  use  v,R  —  u(  and  linear  strain-displacement  equations,  and  assume  4>r  depends  on  displace¬ 
ment  derivatives  through  the  so-called  effective  strain  invariant:  for  loading  in  a  plastic  state.  4»D  *  4>u  4  <J>f. 
Addition  of  (A. 23)  and  4>r  to  potentials  for  the  special  microcracking  model.  (A. 1 7)  and  (A  IR)  respectively, 
enables  one  to  account  for  an  effect  of  microcracking  on  residual  stresses. 

The  results  in  the  Appendix  concerning  the  interpretation  of  O1’  and  (I><  as  mechanical  work.  (A. 14)  and 
(  A.  15),  as  well  as  the  7V  integral  and  energy  release  rate,  are  valid  with  these  extensions  of  the  damage  model 
However,  each  group  in  the  total  potential  <J>  has  to  be  considered  separately  since,  at  any  given  time  and 
location,  some  may  be  for  damaging  processes  and  others  for  constant  damage  processes. 


Resume 


On  developpe  des  methodes  d’analyse  dcs  deformations  et  de  rupture  quasi  statiques  pour  une  classe  de  ntrleux 
visco-elastiques  non  lineaires  et  on  illustre  dcs  applications  types  de  ces  methodes.  Le  choix  de  la  classe  est  dictc 
par  le  comportement  rheologiquc  reel  de  materiaux  monolitiques  ou  composites,  ainsi  quo  par  la  necessile  de 
simplifier  I'approchc  pour  comprcndre  1’effet  des  param£tres  de  base  du  materiau  et  du  continuum  sur  leur 
comportement  vis-a-vis  de  la  croissance  de  la  fissure.  On  discute  cn  premier  lieu  la  pertinence  des  theories  de 
1‘integrale  J  et  du  taux  de  relaxation  d'energie,  dans  le  cas  de  milieux  elastiques  non  lineaires.  On  ctablit  ensuite 
des  equations  visco-elastiques  non  lineaires  ct  on  developpe  les  principes  de  correspondance  qui  permettent  de 
met  (re  en  place  une  relation  simple  entre  les  £tat  meamiques  correspondent  a  des  milieux  elastiques  et  a  des 
milieux  visco-elastiques. 

On  tire  des  principes  unc  base  pour  etendre  la  theorie  de  Pintegralc  J  a  la  croissance  de  fissures  dans  les 
materiaux  visco-elastiques.  L' accent  cst  place  sur  la  prediction  du  travail  mecaniquc  susceptible  d'amorcer  et 
d'entretenir  la  croissance  d’une  fissure  a  I'extremite  de  cclle-ci.  Quelques  exemples  montrent  comment  le 
comportement  a  la  croissance  cst  influence  par  les  proprietes  visco-elastiques  et  par  I'intcgrale  J. 

L/ctude  couvre  le  cas  d'une  fissure  dans  une  couche  mince  presentant  des  proprietes  visco-elastiques 
distinctes  de  celles  du  substrat.  F.n  annexe,  on  presentc  une  theorie  apparemment  originale  sur  les  materiaux 
elastiques  et  visco-elastiques  comportant  des  microstructures  evolutivcs  et  on  indiquc  les  conditions  scion 
icsquelles  la  theorie  de  la  rupture  discutee  dans  le  memoire  cst  applicable. 
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Many  important  rooulto  on  tho  dofomotlon  and  fractura 
of  llaaar  and  nanllaaar  alaatlc  aatarlala  hava  boon  obtalaad 
by  uolna  attain  anargy  potantlala  to  charaetarlta  natarlal 
raaponaa.  Saaldaa  a*  nr  log  aa  tho  baala  for  powerful  nathodo 
of  oaact  and  approximate  atmctural  analyala,  local  and 
global  attain  anargy  functlona  hava  boon  uaod  In  tha  pradlc* 
tlon  of  affoetlvo  conatltutlva  propartlaa  (or  thalr  upper 
and  lover  bouade)  of  nultlphaae  nedla  In  terna  of  propartlaa 
and  geometry  of  tha  phaeea;  Included  ara  atudlaa  of  tha 
affact  of  conatant  danage, In  tha  form  of  email  dlatrlbutad 
cracke  and  volda, on  the  global  etrese-straln  behavior  of 
monolithic  and  composite  naterlala  [e.g..  1,2]. 

Strain  energy  concepts  ara  usad  widely  In  fractura 
nachaalcs  [j].  For  alaatlc  natarlals,  tha  nechaalcal  work 
available  at  a  crack  tip  for  producing  crack  growth  la 
equal  to  tha  decrease  la  potential  energy  (consisting  of 
global  strain  energy  and  the  boundary  work  potential),  and 
this  relationship  has  resulted  In  remarkably  auccaeeful  in¬ 
vestigations  of  fractura  of  rubber  In  its  nonlinear  range 
of  behavior  [*],  ae  well  aa  linear  elastic  naterlala. 

In  this  talk  we  discuss  Inelastic  behavior,  and  con¬ 
sider  Che  question  of  whether  or  not  potantlala  which  are 
enalogoue  to  strain  energy  nay  be  used  to  characterise  de¬ 
formation  and  fracture.  Tha  Ideas  are  Illustrated  using 
date  on  polynsrlc  composite  materials.  First,  elastic 
naterlala,  but  with  timewise  conatant  and  changing  dis¬ 
tributed  "nlcrodsmagc" ,  are  discussed.  Then  viscoelastic 
behavior  and  macrocrack  growth  are  considered  briefly. 

Theory.  The  concept  of  a  potential  for  elastic 
materials  with  damage  nay  be  Introduced  through  the  uniaxial 
stress-strain  curve  In  Fig.  1.  Let  us  suppose  that  a  pre¬ 
viously  undamaged  specimen  Is  strained  nonotonlcally  until 
the  strain  Is  ca<  The  strain  Is  then  reduced,  as  shown  In 
Fig.  1.  Assuming  that  tha  bar  Is  elastic  end  has  constant 
damage  during  the  unloading  period,  with  stress  ou,  and  that 
the  maximum  strain  c,  serves  to  define  the  amount  and  effect 
of  damage,  we  nay  write 

ou  -  fU.c^)  (1) 

Tht  «i»Uw  •train  at  any  t&JM  aquala  tha  currant  strain  on 
tha  loading  cunra,  and  thus 

9 •  f(c,c)  (2) 

Tha  nachanlcal  work  (par  unit  initial  voIim)  during  loading 

to  an  arbitrary  strain  Is 

♦l  •  ♦*•(«)  f  ArL  dt”  -  / f(c\  c')dc'  (3) 

o  n 

where  the  prime  denotes  a  dimmy  variable  of  Integration. 

Tha  net  work  Input  to  tha  ssmpla  at  any  time  during  unload¬ 
ing  la  the  a  haded  area  In  Fig.  1, 

♦u  -  «U(c,c_)  s  ♦t(t_)  ♦  /*  oU  dt' 


*  *t(«m>  *  £  <«'  (*> 

Observe  that  during  loading  and  unloading,  raspactlvsly, 

o1  -  d*L/dc,  oU  •  »«U/»c  (3) 

It  la  convenient  to  lat  ♦  denote  a  quantity  which  aquala 
during  loading  (tg*  «)  and  equals  du  during  unloading 
(c<<a).  Then,  we  nay  write  for  both  loading  and  unloading 
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processes, 

o  •  )«/)■  (g) 

and  tharafora  tha  quantity  ♦  la  a  "work  potential"  which  be¬ 
comes  tha  usual  strain  anargy  density  whan  tha  loading  and 
unloading  curves  ara  Identical.  Obviously,  a  work  potential 
#  can  always  be  constructed,  given  tha  uniaxial  atraaa-atraln 
curves,  Eqa.  (1)  and  (2). 

For  characterisation  of  tha  mult lax lal  atraaa-atraln 
behavior,  or  for  other  raaponaa  functions  which  depend  on 
ears  than  one  Independent  Input,  n  work  potential  does  not 
necessarily  exist.  However,  that  It  can  ba  expected  to 
exist  for  none  realistic  situations  will  ba  dlscuased  nara. 
For  purposes  of  generality,  lat  us  use  as  Independent  In¬ 
puts  tha  generalised  displacements  q,  (]*1,2,...J).  The  re¬ 
sponses  are  the  generalised  forces  K,  which  ara  defined  la 
the  normal  way  by  tha  condition  that,  for  each  j, 

TO  -  QJdq1  (7) 

where  dW  la  the  virtual  work  Input  associated  with  the 
virtual  displacement  <q,.  Suppose,  for  example,  we  let  each 
q.  represent  an  Independent  component  of  a  suitably  defined 
three-dimensional  strain  taaeor  and  let  TO  be  virtual  work 
par  unit  Initial  volume.  Then,  J-6  and  Eq.  (7)  implies  the 
set  Qj  represents  the  components  of  a  stress  tensor  (for 
large  or  small  strains).  In  order  to  characterise  tha  be¬ 
havior  of  laminates,  one  nay  want  to  Identify  the  aat  q, 
with  the  middle  surface  curvatures  and  strains.  In  this 
case,  the  Q.  would  correspond  to  moments  and  ln-plane  forces 
per  unit  length. 

As  in  the  uniaxial  example,  we  assume  that  when  the 
damage  la  constant  the  body  (material  element,  teat  speci¬ 
men.  or  complete  structure)  is  elsstic  in  the  ususl  sense, 
l.e.,  a  strain  energy  function  or  work  potential  Se  exists 
with  the  property  that 

Qj  -  »*c/»qj  (S) 

(gather  chan  using  the  terms  "loading"  and  "unloading"  we 
shall  instead  now  refer  to  "damaging  processes"  and  "con¬ 
stant  damage  procesaes",  since  we  do  not  want  to  imply 
that  the  damage  is  always  constant  when  the  magnitude  of 
one  or  more  loads  or  dlsplacaments  decreases  with  time.) 

The  effect  of  damage  on  Q,  is  assumed  to  be  fully  repre¬ 
sented  by  a  set  of  parameters  Fcn  (n-1,2 . N) .  Following 

the  arguments  in  [3]  it  can  be  shown  chat  a  work  potential 
I®  exists  during* damaging  processes  such  that 

Qj  -  >«D/><ij  (9) 

where  I®  li  a.  function  of  only  the  current  veluee  of  q.,  If 
♦c  ha*  the  font*  ^ 

N 

♦e  *  *  £  '(f.u  «o) 

o  i  n  n  cn 

n-i 

where  *0  la  the  work  potential  without  damage.  The  functlona 
Fn"Fn(q,)  are  to  be  selected  such  that  Fcn”  Fn  during 
damaglngJ processes.  (To  recover  the  uniaxial  case,  Eqs.  (1) 
and  (2),  we  set  N-l  and  cake  F.«q  «C,  F  •  c^) .  1*  should 

be  noted  that  "Internal  variables"  are  hot  ueed  in  thla 
theory.  Of  course,  one  could  think  of  the  set  Fcn  as  "in¬ 
ternal  parameters”  which  are  conatant  during  conatant  d — g- 
procasses  and  which  vary  in  such  a  way  that  Fc_  •  Fn  for 
damaging  procesaes.  For  laminates,  an  Fn  may  be  an  invariant 
of  ply  (or  ply-pair)  strains,  where  the  summation  in  Eq.  (10) 
would  extend  over  all  piles. 

Ceneralliatlone  of  thla  formulation  are  discussed  in 
[3].  For  example,  at  any  given  time  some  of  the  terms  ♦ 

In  Eq.  (10)  nay  be  for  constant  damage  processes  while  others 
are  for  damaging  processes.  Also,  tha  potentials  may  depend 
explicitly  on  time,  and  thus  provide  for  affects  of  sglng  or 
changing  physical  environments.  (The  theory  In  [3],  which 
allows  for  large  deformations,  uses  dlsplacsnent  gradients 
sod  Flols  stressss  instead  of  generalised  dlsplacaments  and 
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fores*.  However,  the  earlier  formulation.  Including  ice  ex¬ 
tension  to  vlscoelascic  behavior,  carrye  over  fully  In  earns 
of  Che  generalized  variables  used  here.)  Thus,  as  In  Eq. 
(6),  we  may  wrica 

Qj  -  3*/3 qj  (11) 

even  If  Che  daeage  parsaecers  In  soma  components  »n  are  con¬ 
stant  and  others  are  varying  in  else.  Also, 


(l.J  -  1.2, ...J) 


except  at  the  points  of  change  froa  one  process  to  another, 
considering  all  *,,;.  The  derlvacivee  in  Eq.  (12)  are.  In 
general,  discontinuous  at  these  transition  points  (cf.  Fig. 

1)  and  thus  Eq.  (12)  does  not  apply  there;  we  assuae,  how¬ 
ever,  chat  the  derivatives  are  continuous  functions  of  dis¬ 
placements  at  all  other  points.  Evidence  of  transition 
points  may  appear  In  experimental  data  as  non-zero  values 
of  this  difference  of  derivatives  (for  1PJ)  over  short  time 
Intervals. 

It  Is  not  necessary  for  the  constant-damage  potential 
to  have  the  form  in  Eq.  (10)  for  the  work  potential  In  Eq. 

(11)  to  exist.  For  example,  a  different  form  for  >c  was 
given  in  [ 5]  which  contains  the  deformation  cheory  of  plas- 
ciclcy  (wlch  elastic  unloading).  We  may  extend  the  present 
model  to  accounc  explicitly  for  stable  micro-  or  aacrocrack 
growth.  (Certain  microcrack  effects  are  at  least  Implicitly 
contained  in  the  *  discussed  so  far.)  An  important  extension 
is  obtained  by  adding  a  term  which  is  equal  to  the  mechani¬ 
cal  energy  required  at  crack  tips  to  extend  any  micro-  or 
macrocracks  not  already  accounted  for  in  9.  For  example,  if 
che  energy  required  per  unit  of  surface  area  A  (projected 
onto  the  local  crack  plane)  were  constant,  with  value  r,  then 
one  would  add  Che  "fracture  potential"  TA  to  *  to  obtain  che 
total  potential.  Without  assuming  che  fracture  energy  T  la 
constant  or  is  che  same  for  all  cracks,  let  us  suppose  a 
fraccure  potential  *f  -  9f(Ay,t)  exists,  where  Ay 
(k*l,2,...K)  represents  che  sec  of  all  (oriented)  crack  sur¬ 
face  areas  needed  to  define  the  effect  of  all  cracks  not 
accounted  for  in  *;  time  t  is  shown  to  indicate  chat  che 
required  crack-tip  work  may  be  affecesd  by  material  aging, 
f-tlgue.  transient  temperature,  etc.  The  coca!  potential 
is  * 


where  «  is  chat  in  Eq.  (ID.  The  mechanical  work  available 
ac  che  edge  of  che  kth  crack  area  Ay  is  -3*/3Ay,  where 
*”*(qj.A.).  whether  or  not  che  damage  parameters  Fca  vary 
U.  P-22z].  Thus,  for  stable  quasi-static  crack  growth  or 
no  growth,  respectively, 

3*{/3Ay  -  -3*/3Ak  or  3*f/3Ay  >  -3*/3Ay  (14) 

The  first  relationship  is  a  set  of  equations  for  finding  all 
of  che  corresponding  A.  as  functions  of  q,.  Using  both  cases 
in  Eq.  (14),  it  is  readily  shown  chat 

Qj  -  3*r/3qj  (15) 

with  cooacanc  crack  areas,  and  wlch  growing  cracks  when  che 
derivatives  3Ay/3qj  exist.  All  Ay  for  mlcrocracks  which 
are  uniquely  determined  as  functions  of  q,  by  che  first  ex¬ 
pression  in  Eq.  (14)  may  be  interpreted  as  damage  parameters 
and  thus  Included  (possibly  as  functions  of  Ay)  in  ths  sec 
Fcni  the  fracture  energy  of  these  Ay  would  not  be  contained 
In  9, .  When  multiple  solutions  for  growing  Ay  exist,  che 
displacement  history  is  needed  to  be  able  to  predict  che 
actually  occurlng  values.  For  this  case  and  for  macrocracks 
one  could  not  Include  Ay  (or  functions  of  A^)  in  the  Fcn. 

Experimental  Studies;  Farris  [6]  showed,  la  effect, 
chat  a  particle-filled  rubber,  wlch  void  growth,  under  axial 
stress  and  confining  pressure  obeyed  Eq.  (12)  (with  q  •  c, 
q2  •  dilatation)  for  constant  axial  strain  rates.  Thus,  at 
least  in  a  limited  study,  the  work  potential  with  damage 
growth  was  established.  Strain-history  effects  due  to  vis¬ 
coelasticity  and  damage  have  been  taken  lnco  accounc  in 
more  recent  work  [7],  We  are  presently  studying  che  valid¬ 
ity  of  the  elastic  and  viscoelastic  theory  with  damage  for 


graphlte/epoxy  tube  and  place  specimens.  In  some  esses 
there  is  a  very  significant  nonlinear  effect  of  axial  strain 
on  shear  behavior,  as  Illustrated  In  Fig.  2.  Specimen  fail¬ 
ure  often  occurs  by  delamlnatlon,  and  this  type  of  macro¬ 
crack  la  being  analyzed  using  Eq.  (14).  Recent  findings 
from  this  investigation  will  be  described. 
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Figure  1.  Uniaxial  scress-ecraln  curve  for  elastic  material 
wlch  increasing  damage  during  loading  (oL)  and 
constant  damage  during  unloading  (ou) . 
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Figure  2.  Shear  stress-strain  curves  with  changing  axial 
strain  c,  from  axial-torsion  tests  of  flat  rec¬ 
tangular  laminates  of  fiber-reinforced  resin. 
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ABSTRACT 

Methods  of  quasi-static  deformation  and  fracture  analysis  have  been  developed  for 
nonlinear  elastic,  viscous,  and  viscoelastic  materials  with  distributed  damage  [l].  The 
crack  growth  theory,  which  uses  a  generalized  J  integral  that  allows  for  viscoelasticity 
and  distributed  microscale  damage,  is  not  much  more  involved  than  that  of  nonlinear 
elasticity  or  special  cases  of  linear  viscoelasticity.  This  simplicity,  compared  to 
what  one  may  expect,  is  a  direct  result  of  the  particular  type  of  constitutive  equations 
and  mechanical  variables  selected  to  characterize  rheological  behavior.  Considering 
elastic  materials  with  distributed  damage,  for  example,  the  constitutive  theory  is  ex¬ 
pressed  in  terms  of  one  strain  energy-like  potential  for  loading  and  another  for  unload¬ 
ing.  The  research  activity  is  presently  in  the  early  stages  of  an  investigation  of  the 
applicability  of  the  theory  to  deformation  and  fracture  of  fiber-reinforced  plastics. 

It  is  anticipated  that  additional  information  on  its  applicability  will  come  from  other 
AFOSR-sponsored  projects  at  Texas  A&M  under  the  direction  of  Professors  Allen,  Bradley, 
Kinra,  and  Weitsman. 

In  this  presentation  we  illustrate  some  features  of  the  theory  for  elastic  compos¬ 
ite  materials  with  damage  and  discuss  current  research  activities.  One  important  ques¬ 
tion  is  concerned  with  whether  or  not  strain  energy-like  potentials  actually  exist;  it 
has  been  addressed  theoretically  in  [1]  with  encouraging  results.  Some  examples  of  real 
nonlinear  material  behavior  which  can  be  characterized  in  this  manner  are  given  in  this 
presentation,  and  then  our  experimental  program  to  study  this  question  for  nonlinear  be¬ 
havior  of  unidirectional  laminates  is  described.  Another  portion  of  the  work  is  con¬ 
cerned  with  application  of  the  theory  to  fracture  characterization  and  analysis, 
assuming  the  requisite  potentials  exist.  In  this  case  the  finite  element  method  is 
being  used  to  predict  crack  initiation  and  growth  in  materials  with  large-scale  dis¬ 
tributed  damage,  first  for  initially  isotropic  and  homogeneous  media  and  then  for  com¬ 
posites  (i.e.,  delamination  initiation  and  growth).  A  few  years  ago  we  employed  linear 
elastic  fracture  mechanics  in  an  investigation  of  the  fracture  behavior  of  a  randomly 
oriented  glass  fiber  reinforced  plastic,  SMC-R50  [2].  The  data  on  this  nonlinear  mate¬ 
rial  are  reinterpreted  here  using  J  integral  theory  in  order  to  further  illustrate  its 
use.  Application  of  J  integral  theory  to  delamination  growth  when  large-scale  distrib¬ 
uted  damage  exists  is  under  study  on  Professor  Bradley's  project. 
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CONSTITUTIVE  THEORY  USING  POTENTIALS  FOR  MATERIALS  WITH  DAMAGE 
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Uniaxial  Loading 
aL  =  f(e.e)  =  ^ 


Multiaxial  Loading 

aL.  =  3 <t> L/ 3 e  .  . 

13  13 


Unloading 


oU  =  f(e»W  =  ^- 


Unloading 
°ij  =  3*U/3eij 


FEATURES 

•  CONVENIENT  FORMULATION  FOR  FRACTURE  APPLICATIONS 

•  DAMAGE  PARAMETERS  (E.G.  £„„)  ARE  EXPRESSED  IN  TERMS  OF  STRAIN  HISTORY  AND  MAY  BE 

ITiaX 

RELATED  TO  PHYSICS  OF  DAMAGE  PROCESS 

•  POTENTIALS  ANALOGOUS  TO  STRAIN  ENERGY  DENSITY  ARE  USED  FOR  LOADING  (4>L)  AND  UNLOADING 

(*°) 

•  ALLOWS  FOR  TEMPERATURE  AND  MOISTURE  INDUCED  STRESSES 

•  VISCOELASTIC  EFFECTS  ARE  INTRODUCED  BY  USING  "PSEUDO  DISPLACEMENTS"  IN  PLACE  OF 
DISPLACEMENTS 

APPLICATIONS 

•  EXISTENCE  OF  POTENTIALS  SHOWN  IN  THE  SPECIAL  CASES  OF  PARTICLE-REINFORCED  RUBBER, 
ELASTO-PLASTIC  BEHAVIOR  OF  METALS  AND  SECONDARY  AND  TERTIARY  CREEP  OF  METALS 


•  USE  FOR  FIBER-REINFORCED  PLASTICS  IS  UNDER  STUDY 
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CONSTITUTIVE  EQUATION  EXAMPLE 


A  SPECIAL  PARTICULATE  COMPOSITE  WITH  MICROCRACKING 
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<fL  =J  -  g(W')JdW' 

fWM 

*U  =  [l  ~  9<WM)][W  -  «M]  +  J  [l  -  g(W)]c 

where  Vi  =  STRAIN  ENERGY  DENSITY  WITHOUT  DAMAGE  (g  =  0) 
WM  =  MAXIMUM  W 

•  NOTE  THAT  «L  =  $U  WHEN  W  =  W„ 

M 

•  STRESSES  ARE 
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■GIVES  REDUCED  STIFFNESS 
WHEN  0  <  g  <  1 
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APPROACH  FOP  UNIDIRECTIONAL  PLY  CHARACTERIZATION 
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CONSIDER  A  UNI-AXIALLY  LOADED,  OFF-AXIS,  UNIDIRECTIONAL  COMPOSITE 
TENSILE  SPECIMEN.  THE  STRESSES  IN  THE  PRINCIPAL  MATERIAL  DIRECTIONS 
ARE  WRITTEN  IN  TERMS  OF  THE  STRAINS  IN  THOSE  DIRECTIONS  AS  (NEGLECT¬ 
ING  EFFECTS  OF  END  CONSTRAINT) : 


Oj  -  o1(c1  ,  c2  ,  Y12> 


°2  *  °2<£1  •  £2  *  Y12>‘ 


IF  A  POTENTIAL  ♦  IS  TO  EXIST,  SUCH  THAT 


,  o j  -  a*/ac2  ,  t12  - 


IT  IS  NECESSARY  THAT,  FOR  EXAMPLE, 
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•  PERFORM  TESTS  ON  AS4/3S02  MATERIAL  USING  OFF-AXIS  SPFCIMENS. 
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0  *  FIBER  ANGLE 

BY  VARYING  0,  THE 
NECESSARY  DATA  CAN  BE 
DEVELOPED  TO  EVALUATE 
THE  CPOSS- DERIVATIVES 
AND  DETERMINE  THE  EXIST¬ 
ENCE  OF 


THESE  DATA  ALSO  CAN  BE  USED  TO  DEVELOP  THF  POTENTIAL  t  ITSELF. 
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FRACTURE  EXAMPLE 


Measured  by  o 
45,000  N-copccity 
load  cell 
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Figure  3.  Double  edge  notch  specimen  with  measured 
parameters  indicated.  (SMC-R50) 
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Nominal  thickness 
of  all  specimens 
was  2.54  mm 


Figure  4.  Double  edge  notch  specimen  geometry.  (SMC-R50) 


Figure  5.  Stress-strain  curves  based  on  load-displacement  data 

for  25.4  nm-wide  unnotched  tensile  specimens.  (SMC-P50) 


ANALYZE  THE  INITIATION  AND  PROPAGATION  OF  A  CRACK  IN  A  MATERIAL  WITH 
DAMAGE  THROUGH  THE  USE  OF  A  FAILURE  ZONE  CRACK  TIP  MODEL  AND  THE 
GENERALIZED  J  INTEGRAL. 


INITIALLY,  THE  ISOTROPIC  MATERIAL  MODEL  OF  J2  DEFORMATION  THEORY  OF 
PLASTICITY,  MODIFIED  BY  ELASTIC  UNLOADING,  IS  BEING  USED*. 

°lj  *  2Wcij  '  for  T«  ‘  (t»’n»x 

°ij  *  aii  =  2G<cij  '  ‘ij1  '  for  T*  *  (Vmax 


WHERE:  0^.  •  o. }  -  j 

xs  -  i^orcr^2 


EXTENSION  TO  ORTHOTROPIC  COMPOSITES  WITH  DAMAGE  TO  FOLLOW 


•  FAILURE  ZONE  MODEL  FOR  THE  CRACK  TIP: 


failure 


r  -  fracture 
energy 


•  THE  J  INTECRAL  IS  USED  TO  OBTAIN  INFORMATION  ABOUT  WORK  INPUT  TO  THE  FAILURE 
ZONE  FROM  REMOTE  FIELD  QUANTITIES: 


Figure  8.  Deflections  in  a  cantilever  supported  by  an  elastic-plastic  foundation  in 
which  the  tractions  vanish  above  a  critical  displacement. 


x  fin.) 


Figure  9.  Crack  opening  displacement  profile  for  thick,  nonlinear  bend  specimen. 
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